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Abstract 

The method of analytic continuation is used to find exact integral 
equations for a selection of finite-volume energy levels for the non¬ 
unitary minimal models M. 2 , 2 N+z perturbed by their <pis operators. 
The N =2 case is studied in particular detail. Along the way, we find 
a number of general results which should be relevant to the study of 
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1 Introduction 


The finite-volume spectrum of a quantum field theory encodes a great deal 
of information, of interest even if the ultimate objective is to understand the 
behaviour of the theory in a world of infinite extent SI- In the unusually 
tractable examples provided by integrable theories in 1+1 dimensions, it 
has been known for some time that an effective way to obtain the lowest 
level of this spectrum (the ground-state energy) is provided by the integral 
equations of the thermodynamic Bethe ansatz, or TBA [3], More recently, 
modifications to these equations have been found which extend the treat¬ 
ment to other (‘excited’) energy levels, for the rather simple but nonetheless 
nontrivial case of the scaling Lee-Yang model, or SLYM Si- (An alter¬ 
native approach, based on the Destri - de Vega equation, was initiated in 
refs. [||], but its relation with the TBA remains to be understood.) In this 
paper we explore how the methods of J5| can be generalised to a series of 
models, namely the perturbations of the minimal models -Ad 2 , 2 A +3 by their 
operators. Apart from IV=1, which is just the scaling Lee-Yang model, 
these theories all have more than one mass in their spectrum, and their 
excited-state TBA equations exhibit a number of new phenomena. The 
bulk of the paper concerns the N =2 case, for which we have some detailed 
numerical results. For higher values of N, our results are less complete, 
though we can be confident that we have captured the exact behaviour of 
the one-particle levels in the infrared region. 

The paper begins with a brief review of the relevant perturbed con¬ 
formal field theories, and a discussion of some general properties of TBA 
equations that will be important later. Ref. || can be consulted for some 
further background on the idea of analytic continuation as applied to TBA 
systems. In section 3 the method is used to give a reasonably complete 
analysis of the first few levels for N= 2, whilst section 4 describes some 
features which emerge when the same techniques are applied to the other 
theories. None of these models is of overwhelming physical importance, and 
our main aim has been to use them as a testing-ground for the whole an¬ 
alytic continuation approach to generalised TBA equations. With this in 
mind, the concluding section 5 includes a summary of those features which 
have emerged during our investigations that should have wider applicability. 
There are two appendices: the first recalls some field-theoretic predictions 
for infrared asymptotics, and the second compares the numerical solutions 
of our equations with results obtained from the truncated conformal space 
approach. 
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2 The T/v-related perturbed conformal field theo¬ 
ries 


The minimal models M. 2 ,2N+3 have JV+1 primary fields , occupying a 
single row of the Kac table. The central charge is c = — 2N(6N+5)/(2N+3), 
and the scaling dimensions of the primary fields are 

(s—l)(2iV+2—s) 


di, s = -- 


s = 1...N+1 . 


2JV+3 

Perturbing by ipis results in a factorised scattering theory sometimes de- 
noted by T/v, or else 0,B- This has N self-conjugate particles with 
masses M a : 


sin( na/h) 

M a = ;,/ Mi 


a = 1...N 


sin(7r /h) 

where h = 21V+1, and two-particle S-matrix elements S ab (0 ): 


a-\-b— 1 


Sab= n (a, 6 = 1 ... A) ; 


( 2 . 1 ) 


|a —6| + 1 
step 2 


{x} = (x — 1 ) (x + 1 ) 


(.)(*)==» 


17VX > 
2/t > 


sinh(| — 


Z7ra; > 

2h > 


The exact relation between the <p 13 coupling A and the mass Mi of the 
lightest particle was found by Al.B.Zamolodchikov ['9], Setting g(x) = 

v/r(x)/r(l-s), it is 


4sin£ 

Mi (A) = -A-- 


7 r 


’(s+h) 


7T 


(h 2)(h 4) /h- 4 \ 

(/r+2) 2 9 \h+2)9{ h+2 ) 


h +2 
4/i 


h+2 

(—A) 4h 


( 2 . 2 ) 

In the normalisations implicit in [9j, a sensible scattering theory with real 
masses for the particles is obtained when the coupling A is real and negative. 
The only exception is N= 1, where the imaginary values of some conformal 
structure constants mean that a real perturbative expansion is obtained if 
A is purely imaginary, and ImA must be positive for real masses in the 
infrared 


10]. 


Since the S-matrix (2.1) is diagonal, the TBA equations for the ground 
state energy on a circle of circumference R are §0 : 

N 

£ a (d) = m a r cosh 6 - ^ cj> ab *L b (9) , (2.3) 


6=1 
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with 


r = M 1 R, m a = ^, 4>ab(9) = log S ab (9), L a (0) = log (l+e £a(0) ) , 

and f*g(9) denoting the convolution -Y I^° 00 d9' f(9--Q')g{0') . Solving these 
equations yields N pseudoenergies e a (9 ); the ground state energy is then 

E 0 (X, R) = £Wik(A, R) - j ^c(r) ; 

bR 

^buik(A, R) = - M )[ X} * r , c(r) = 4E [°° d9r cosh 0L *( e ) ' ( 2 ‘ 4 ) 

8 sm(27r/ n) ^ ~ l J-oo 

It will sometimes be convenient to work with the scaling functions Iq(r), 
related to the energy levels E % by Ei(X,R) = ^ Fi(M\R ). 

The functions Y a (9) = exp(e a (0)) provide an r-dependent set of solutions 
to a set of functional relations called a Y-system |l2|, |l3| 

N yjp j 

Y a (9 - f)Y a (9 + f) = I] (! + W))^ , (2.5) 

6=1 

It 1 

where l l ab Ni = 5 a}b _i + 6 a< b+i + 5 a ,Ndb,N is the incidence matrix of the T/y 
graph. A (none-too-obvious) consequence of these equations is that the 
functions Y a (9) are wr(h+2) //i-periodic: 

Y a {9 + ^d) = Y a {9) . (2.6) 

Subtracting the right-hand side of ( |2.5| ) from the left-hand side yields an 
analytic function of r and 9 which is identically zero. This function therefore 
remains zero during any process of analytic continuation, and so the Y- 
system and all of its consequences hold equally for all the continuations of 
the basic TBA to be encountered below. With this in mind, it is worth 
pausing to record a couple of other general properties before proceeding to 
specific examples. 

The first concerns the behaviour of the Y a as functions of r and 9 com¬ 
bined. If we define a new pair of variables as follows: 

a. |_ = [re u J ; a_ = [re (2.7) 

then we claim that each Y a (a + , a_) can be expanded as a power series about 
a+ = a_ = 0, with a finite domain V C ® x (E of convergence. This is a 
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conjecture, but note that it corresponds to an expansion about the central 
region in the ‘kink limit’ of the TBA, where we would not expect to find 
any singularities, all pseudoenergies tending to constants there. 

The second general property follows from the first. Consider simultane¬ 
ous changes of r and 9 which leave a + and a_ invariant. These take the 
form r —> e ia r ; 9 — > 9+i(3 with a = (p+q)n(h+2)/2h, (3 = (p—q)7r(h+2)/2h, 
and (p. q) a pair of integers. The (—1,0) case of this transformation shows 
that 

y o ( e ~* w ( fc+2 )/ 2/ V, = Y a (r , 9) ((a+, a_) € V). (2.8) 

This result is only guaranteed while |rj and \9\ are small enough that (a+, a_) 
lies in V. However as functions of 9 alone, the Y a are expected to be entire 
(in the conformal case this can be proved in some situations j|). If |r| 
is sufficiently small, then (a+,a_) will lie in V for, say, real 9 in the unit 
interval; the identity (|2.8| ) holds for these values of 9 , and can then be ana¬ 
lytically continued to all 9 while r remains fixed. Using the real analyticity 
of the Y's we can then deduce that for these sufficiently small values of |r| 
the functions Y a (r,9 ) have the symmetry 

Y a (r,9) = Y a *(r,6) (2.9) 

where 

r = e Mr(A+ 2)/ 2 V , 9 = in^+9* (2.10) 

and the asterisks denote complex conjugation. For larger values of |r|, the 
point (a+,a_) will be outside V for all 9, and the argument breaks down. 
Indeed, the functions Y a (r, 9) acquire a multi-valued structure reflecting the 
multivalued nature of the energy levels themselves. However, since all of 
these functions can be found by analytic continuation from the small-|r| 
region, a memory of the earlier relation remains in an involution on the 
set of all T’s: given a function Y a (r,9 ), the exponential of the solution to 
some (perhaps generalised) TBA equation, another function Y a (r,9) can be 
defined by 

Y a (r,9) = Y a *(r,9) . (2.11) 

The logarithm of this function, £ a (6) = logy a (r, 9) , will again solve a (pos¬ 
sibly different) TBA equation, with the radius r replaced by the reflected 
radius r. The observation turns out to be rather useful, allowing solutions 
of the generalised TBA equations to be obtained in some situations where a 
direct numerical attack fails. So as not to have to rely for it on the power- 
series expansion, we also checked its validity, for N=2, 3 and 4, in various 
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situations where the equations related to both Y a and Y a could be solved 
numerically. 


3 Excited states for Ml 2.7 + ^13 


In this section we restrict our attention to the case N=2, h=5. The approach 
proposed in j5) requires the continuation of r along some path in the complex 
plane, in the hope of inducing a monodromy in c(r ) which will convert the 
ground state into an excited state. To decide which paths to take, we used 
the truncated conformal space approach (TCSA) [10], and in particular the 
program of ref. |l4|], to map the first few sheets of the Riemann surface 
associated with the ground state scaling function Tb(r). The results are 
shown schematically in figure [l| 



Figure 1: The first five sheets of the Riemann surface of the scaling function 
F(r), showing on the left the connectivity of the various sheets, and on the 
right the positions of the square-root branch points over the complex r-plane. 

Figures 2 and 3 may also help to visualise how the sheets fit together. 
Figure 2 plots various F(r ) for real values of r between 0 and 12, labelling 
the lines both by their ultraviolet limits near r=0, and by their asymptotic 
mass gaps at large r. Figure 3 shows the real parts of these same F(r ) along 
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Figure 2: TCSA data for real r between 0 and 12, with the lines labelled 
both by the conformal fields found in the ultraviolet, and by the mass gaps 
exhibited as r —> oo. 


F(pe 7 ” i ' 20 ) 



Figure 3: TCSA data for complex r on the positive-A line r = pe 7 ™/ 20 , 
0 < p < 12. Three branch points can be seen; those marked A and C also 
appear on figure 1. 
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another line of interest, namely r = pe 7m / 20 , p E 1R + . Via equation (2.2) and 
the relation r = this line also corresponds to real R and A, but now 

with A positive. Consideration of the power series expansions in A provided 
by perturbed conformal field theory shows that the scaling functions must 
be real here too on some initial segment of the positive-A line. Beyond the 
radius of convergence of the perturbative expansions singularities are found, 
after which the F(r) can pick up imaginary parts. Examples of singularities 
on the positive-A line are the branch points marked A and C in figures 1 
and 3. 

As noted in ||, another consequence of the perturbative series is that all 
other branch points occur in pairs r, r, related by a complex conjugation of 
A, which reflects the complex r-plane in the positive-A line. This parallels 
the involution of the Y’s discussed at the end of the last section. For N= 2, 
the relationship is 

r = e 7m / 10 r * . (3.1) 


Examples in figure 1 are the pairs B, B] D, D ; and E, E. 

Using the branch points, it is possible to move between the energy levels. 
A path encircling A will link the ground state with the first excited state, 
one encircling A and then B will link the ground state with the second ex¬ 
cited state, and so on. During this continuation, singularities in one or other 
of L\ and L 2 may cross the real axis. Whenever this happens, all integra¬ 
tion contours in the TBA equations must be distorted to run around these 
singularities. When the contours are returned to the real axis, residue terms 
are picked up, resulting in modified equations in which the positions of some 
singularities make an explicit appearance. We will call these singularities 
‘active’, and the others ‘inactive’. 

This was described for the case of the SLYM in [S]. The principles are 
the same here, but the practice is considerably more complicated, and the 
rest of this section will be devoted to a detailed description of how it goes for 
the cases of zero, one and two pairs of active singularities (note, the 0^—6 
symmetry of the initial TBA means that singularities always either appear 
or disappear from the equations in pairs). It will sometimes be convenient 
to use the functions z a (9), a = 1,2, where 


z a (0) = 1 + e~ £a ^ = 1 + Y a (d)~ 1 . 


Singularities of L a (0) = log z a (9) are found at zeroes and poles of z a , where 
Y a takes the values —1 and 0 respectively. It will be important later that 
these singularities never come singly: the T 2 Y-system (2.5) ties together 
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the values of the Y a {6) at locations separated by integer multiples of m/5. 
We will label the zeroes of Y] (6) as Ypp, j = 0... oo. With each zero, the 
following other special values are forced: 


k 

-1 0 1 

Yi((3p + 
y 2 (pP + 

y 0 l/y ... 

-1 0 -1 


Here y is free, although it should not be zero since, as mentioned earlier, Y\ 
and Y 2 are expected to be entire functions of 9. The zeroes of Y 2 (9) occur 
either at the points ±/?2 already captured in the above table, or else at a 
further set of points Ypp ' 1 , which entail a different singularity pattern: 


k 

-1 0 1 

Y^pp + i&) 
Y 2 {PP + lJ f) 

— 1 a —1 

y 0 (1 +a)/y ... 


This time both y and a are nonzero (the latter condition preventing singu¬ 
larities in some of the omitted terms in the table). 

As for the other type of singularity in L a (9 ), corresponding to a zero of 
z a (9) with Y a (6) = —1, it is not hard to check that Z2(0) = 0 forces 
with pp equal to either 9—m/5 or 9+m/5, while z\(6) = 0 leads to 


this time with Pp equal to either 9—m/5 or Q+m/5. This is the reason 
for the at-first-sight perverse labelling: for the analytic continuation of the 
TBA equations, it is the zeroes of the z a (9), located at PP Ym/5, that turn 
out to play the crucial role. 

To unravel the singularity structure of any particular solution, it often 
helps to reconstruct its ‘history’ through the analytic continuation, remem¬ 
bering the 9 —► —9 symmetry of the equations. Also useful are the facts that 
whenever r is real, the singularity pattern must additionally be symmetrical 
under the standard conjugation 9 —> 9*, and that whenever r lies on the 
positive-A line, so that r = r in ( |2.10 ), the pattern is either symmetrical 
under the shifted conjugation 9 —* 9, or else is mapped by this operation 
onto the pattern for another solution Y a . 

For later use, we will record the situation for the solutions to the basic 
TBA equation ( |2.3| ), when r is real. A typical pattern, in fact that found at 
r= 1, is shown in figure 0. The points YfdP all have imaginary parts equal 
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Re 

0 12 3 4 



Re 

0 12 3 4 



Figure 4: Contour plots of \z a {9)\/{\+\z a {9)\) in the complex 9 plane, for 
the basic TBA equation with r=l. Concentric patterns of contours are signs 
of singularities in L a {6). Those corresponding to zeroes of Y a (6) are marked 
by dots • ; the others are zeroes of z a (9). The scales on the axes are in units 
of 7t/5. 


to 77t/10 (the same, by the periodicity ( |2.6|) , as — T7t/ 10), and so all of the 
zeroes of Y\ and Y -2 lie on the line Irn 9 = ±77t/ 10. The zeroes of z\ and Z 2 , 
as follows from (3.3) and (3.2) respectively, lie on the lines Irn 9 = ±7r/2. 
For definiteness, we will choose the /3’s to have negative real parts, and to 
be ordered from left to right: 


... < Re < Re ^” 1) < ... < Re (5^ < 0 . (3.4) 

As r -4 0, solutions to TBA equations of this sort split into pairs of ‘kink 
systems’. One starts near — log(l/r) and runs leftwards, the other starts 
near +log(l/r) and runs rightwards, and in the region between the two the 
pseudoenergies are approximately constant. With the labelling choice just 
adopted, the zeroes of Y\ and Yz in the left-hand system are at {/?2 } and 
{/? 2 *\ Pi')} respectively, and those of the right-hand system are at {—(3 
and {—/? 2 *\ — Pi }• Thus the singularities visible in figure |4| all belong to 
the developing right-hand kink system, and move further to the right as r 
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decreases. 

In the deep ultraviolet, the mutual interaction of the two kink systems 
becomes vanishingly small, and they separately solve kink forms of the TBA 
equations, obtained by replacing every occurrence of r cosh 9 either by \re~ e 
(for the left-hand kink system) or by \re e (for the right-hand system). In 
the kink forms of the equations, the dependence on r is trivial: any further 
changes can be absorbed by a compensating shift in 8. This tells us some¬ 
thing new about the singularity locations in this limit: the singularities in 
the left-hand set are at fixed distances from the point 9 = — log(l/r), while 
those in the right-hand set are oppositely placed, at fixed distances from the 
point 6 = log(l/r). To make this more precise, we have 


lim 

r —>0 L' 


(3a ] + log(l/r) 


_ r (j) 

L a 


(3.5) 


with { Ca a set of constants satisfying Imca * = ±77 t/ 10. At large j, the 
real parts of the Ca are given by 

Re 4 j) -log (^( 2 j + 1)tt) ( j > 1 ) . 

For smaller j these real parts are somewhat distorted, though the counting 
(starting at j = 0) is correctly given by this formula. 

During analytic continuation, the singularities will move around, and 
for most of the time we will use the same labels for the singularities of the 
continued solutions. To avoid ambiguity we must remember the particular 
path in the complex r-plane chosen to reach a given solution, since there 
are closed paths on the Riemann surface around which pairs of singularities 
swap over. Some examples will be seen later. 


3.1 The basic TBA in the complex plane 


This case provides the starting-point for the investigation. If r is continued 
up from the real axis onto the positive-A line, then a function F(pe 7nl ^ 20 ) is 
obtained, from the basic TBA equations ( [2.3| ) and (2.4), which is real out 
to p = po = 2.39342(4). This point matches the branch point A found from 
the TCSA. On this segment Y\ and Yi are self-conjugate in the shifted sense 
of equations ( 2.1C| ) and ( |2.11|) : 

Y a (pe^ 0 ,e) = Y:(pe^ 0 ,d), 


0 < P < Po 


(3.6) 


This means that the singularities in L\ and To ai 'e arranged symmetrically 
with respect to the shifted conjugation 9^8 = 7m/10+9* . Here this 
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happens in the simplest possible way, with the zeroes of Y\ and I 2 lying on 
the fixed lines of the mapping, Irn 9 = ±77t/ 20, and the zeroes of the z’s 
symmetrically arranged on the lines Im 9 = ±37r/20 and Im 9 = ±ll7r/20. 
All zeroes in the left half-plane, with negative real parts, lie on the lower 
triplet of lines, and all to the right lie on the upper triplet. Figure || shows 
the situation at \r\ = p = 1. Given the symmetry (3J5), the zeroes must 


Re 

0 12 3 4 



Re 

0 12 3 4 



Figure 5: Contour plots of \z a {0)\/{l+\z a (0)\) in the complex 6 plane, for 
the ground-state solution to the basic TBA equation with r=e 7 W 20 , p=\. 
Labelling as in figure 4. 


lie exactly on the lines just described, if they are not to be doubled up. 
A similar phenomenon was observed in || for the SLYM, although at the 
time it was not understood in terms of the shifted conjugation symmetry, 
and so was only verified numerically. Such a check was also carried out 
here, confirming expectations to a precision of 15 digits along the whole 
segment 0 < p < po, and then showing the singularities moving away with 
a square-root type behaviour as p passed po and the solutions ceased to be 
self-conjugate. 

In the ultraviolet limit, p —■» 0, the splitting into kink systems occurs just 
as for real values of r. Using the result (3.5) for the asymptotic singularity 
positions, equally valid for these complex values of r, it is easily seen how the 
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singularities move during the continuation from the real axis to the positive- 
A line. For example, the zeroes of Y\ and I 2 in the left-hand system move up 
from the line Im0 = —77r/10 to sit on the line Imf? = —77 t/ 20, while those 
in the right-hand set move down from Im0 = +77 t/ 10 to Im0 = +77 t/ 20. 
(If we had instead continued r down into the lower half-plane, these motions 
would have been reversed.) Comparing figures [| and j| should clarify the 
relationship between the patterns before and after continuation up onto the 
positive-A line. 

Of the zeroes of the z a ( 9 ), those at ±(/3a' > +111 / 5 ) now lie closest to the 
real axis, and will be the most important later. Hence we define 

8a ) = Pa j) + W5 • 


} ( 0 ) 


This was for the ground-state solution to the basic TBA equation. But 
there is another solution to this equation on the segment of the positive-A 
line between p = 0 and p = po- This can be obtained by continuing the 
solution already discussed anticlockwise around the point r = poe 7m / 20 . In 
the process, no singularities in L\{9) or L2(9) cross the real axis and so no 
modifications to the TBA equation ( [2.3| ) are required. However the pseu¬ 
doenergies undergo a nontrivial monodromy, and on evaluating the integrals 
in (2.4), a different answer is found which matches perfectly with the first 
excited state on the segment 0 < p < po, which is the line labelled by ip ±2 
in figure ||. 

During the anticlockwise continuation of the ground-state solution, 9^ 
the rightmost zero of Z2{9) in the left-hand kink system, migrates across to 
join the right-hand system, while the matching zero in the right-hand system 
makes the opposite journey. The final pattern again respects the 9^9 
symmetry, but now with a different distribution for the imaginary parts of 
the singularities between the left and right kink systems. The change in 
the pattern of singularities can be appreciated on comparing figure |^, which 
shows the new excited solution on the positive-A line at p= 1, with figure ||, 
the ground-state solution at the same value of p. 

The movement of singularities as the point p = po is approached along 
the positive-A line is shown in figure [7. which plots the real parts of 9^ 
and 02 °'* against log(l/p), both for the ground state and for the first excited 
state. The slopes tend to ±1 as p —* 0, the slopes of —1 being characteristic 
of members of the left-hand kink system, and the slope of +1 acquired by 
the ‘excited’ zero of Z 2 (9) showing that it has joined the right-hand kink 
system. 
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Figure 6: Contour plots of \z a {0)\/{\+\z a {9)\) in the complex 9 plane, 
for the excited solution to the basic TBA equation with r=e 77r */ 20 , p=l. 
Labelling as in figure 4. 


Equation (T5) must be modified to take account of the defection of the 
two singularities, one from left to right and the other from right to left. 
To prepare for more general situations, which may involve further pairs of 
defectors, we replace (|3.5D by 


lim 

r \—>0 


u a ) Pa\ r ) +l0g(l/r) -C<f } 


= 0 , 


(3.7) 


fj) . ... (j) . . 

where fa is +1 if the singularity at (3a is still part of the left-hand kink 
system, but switches to —1 if it has joined the right-hand system. The values 

of the constants Ca are characteristic of each particular kink solution, while 
(i) 

the Va preserve some memory of how the solution has been obtained by 
continuation from the ground state. For example, it appears that only for 
the ground state do all of the Ca ^ have imaginary parts equal to ±77 t/ 10. 
For the excited state now under discussion, line® has become equal to zero, 
and i4°^ is equal to —1. 

It is instructive to see how the analytic treatment of the ultraviolet limit 
must be modified to cope with this excited solution to the basic TBA. The 
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Figure 7: Real parts of 6\ and 9 2 for r on the positive-A line, plotted 


against log(l/p). The lines marked and 9 2 u> are for the ground-state 


3 ( 0 ) 


solution, and those marked and 9 are for the excited state. 


(o) 


usual manipulations of the kink forms of equations (2.3) and ( [2.4| ) give 

2 


6 


(0 ) = -J2 L +(Ca) 

7T * —, 


a= 1 


where, borrowing some notation from of refs. [[L5|] and (l 

L og Ca ( 1+x ) L °g C a ( x ) 


L + (C„)=l [ 

2 JC a 


l+x 


dx 


(3.8) 


(3.9) 


The C a are certain contours in the complex plane, and Log Ca (z) stands for a 
continuous branch of log(z) along the contour C a ; the starting-point of the 
branch will be determined shortly. The modification to the ground-state 
value of c(0) resides in the form of these contours and branches, and care 
must be taken to identify them correctly. The variable change leading to 

is given by x{9)= exp(—ejf nk (0)), 


(3.8) tells us that the contour C a in 
with 6 real and running from +oo (at the start of the contour) to — oo (at 
the end), and e^ mk (0) the (excited) kink solution to the TBA equation. As 
9 —> oo, all pseudoenergies diverge and hence the contours C a always start 
at the origin. In the opposite direction, the kink pseudoenergies tend to 
finite constants and as a result the contours terminate at points x a = 1 /Y a , 
with the constants {Y \, Y 2 ) providing a stationary (^-independent) solution 


to the N =2 case of the Y-system (2.5): 


Y? = (1 + Y 2 ), 
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(i + Fi)(i + y 2 ). 


(3.10) 


Y‘> = 


There are four solutions to this equation: one, (Yi, Y 2 ) = (—1,0), appears 
not to be relevant (it would anyway require e 2 = — 00 ), whilst the others 
can be written as (Yi, Y 2 ) = (y^. y ^) with s = 1, 3,5 , and 


^ sin(sa7r/7) sin(s(a+2)7r/7) 
2 la 


sin(s7r/7) 2 


(3.11) 


These numbers feature in certain sum rules for dilogarithm functions. When 
the dust has settled, it will be most convenient to have these written in terms 
of the Rogers dilogarithm C(y), defined for real y in the unit interval by 




(3.12) 


(In cases when C a runs from 0 to x along the real axis, with all logarithms 
in ( |Q| ) on their principal branches, L + (C a ) = £(x/(1+x )).) The sum rules 
that we shall need are the following: 


C\ 


Vi+!/! 1 ) ) +£ (i+ 1 9 ? ) ) 

~c(l+y^j -£^ 1 + j 4 3) ^ 


4 7T 2 

7~6~ 

6 7 r 2 

2^ 

7 6 


(3.13) 

(3.14) 

(3.15) 


Note that all of the arguments have been adjusted so as to lie in the interval 
[0,1] . There do exist prescriptions to extend the Rogers dilogarithm beyond 
this range, but at this stage it is better for us to keep our options open, since 
the relevant continuations will anyway be determined by the TBA equations. 

For the ground state there is an additional physical requirement that the 
Y’s be positive and the logarithms real, and this selects the s=l solution 
and, via (3.13), the expected result c(0) = 4/7 is recovered. For excited 


states there are no such constraints, and for the state under discussion here, 
s=3. This can be justified either by an appeal to numerical results at small 
|r|, or else by noticing that, as a result of the shifted conjugation symmetry, 
the functions Y o kmk (0) = exp(e kink (0)) are real along the line Im 0 = 2l7r/20, 
and that they each have a single zero on this line, at /3 2 °^. Since as 6 —y 00 
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both are positive (tending to +oo), they must both be negative at 6 = —oo, 
and this picks out the s =3 solution. 

The contours C\ and C 2 join the origin to the points x\ and x 2 respec¬ 
tively, avoiding the singularity at x = — 1 (note that both x\ and X 2 are less 
than —1 for this solution). Examining numerically the small-|r| behaviour 
of the full TBA equation, we found that the correct path to take loops under 
x = — 1 for Ci, and over x = — 1 for C 2 . 

Finally, we should settle the branches of the logarithms. The branches 
of the L a {6) in the TBA are such that L a (+oo) = 0 (otherwise, the integral 
(2.4) would make no sense). This translates into the branch of Log Ca (l + x) 
being such that this function is zero at x = ( 1 . From this and the way that 
the two contours go round the point x = — 1 follows the behaviour at the 
other ends of the contours: 


Log Cl (l+xi) = log 11 T X \| Z7r: 
Log C2 (l+x 2 ) = log |1 T x 2 \ + nr. 


(3.16) 


(Here and below, log x stands for the single-valued branch of the logarithm, 
with Im logx E (— n, 7t] .) The branch of \jOg Ca x can now be determined 
from the stationary form of the kink TBA equation: 


_kink ( 


2 

E- 

6=1 


er~ c (-oo) = v N ab Lf mk (-oo) 


where 


Nab ~ 2vr ./-oo 


<t>ab(0)d0 = 



(3.17) 


(3.18) 


Taking the imaginary part of ( |3.17| ), and using the relations Log Ca (x a ) = 




kink 


(- 00 ) and Log Ca (l+x a ) = L“ nk (-oo), gives 

Log Cl (®i) = log |cci| - ztt; 
Log C2 (x 2 ) = log|x 2 |-z7r. 


(3.19) 


Now c(0) can be calculated, using the strategy advocated in ref. (lj. First 
deform C\ and C 2 to run along the negative real x axis, save for a small semi¬ 
circle below —1 for C\, and a small semicircle above —1 for C 2 . The contours 
can then be split into segments and the integrands reexpressed in terms of 
the single-valued logarithm logx, always with positive real arguments. For 
L + (Ci) this goes as follows. We have 


2L+(Ci) = A-B 
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with 


A = 


Log Cl (l+x) 


dx 


I Ci x 

/' -1 log(l+x) 


/ 0 X 

f 1 log(l— x) 


dx + r Lofc - (1+x) & 

J -1 X 


/o 


X 

1 log(l— x) 

X 


dx + 


dx + 


f xi log(— 1 —x) 


dx — 


fX 1 


i7T 


dx 


/-I x 
f-xi - 1 log X 


/o 


x +1 


J- 1 X 

dx — ITT log( —Xl) 


and 


B = 


C X1 Logci ( x ) 
0 l+x 

r 11 log(—x) 

0 l+x 

/'* logx 

J o 1 x' 


dx 


dx — 


27T 


dx 


/o 


dx — 


l+x 

r-m-i 


/o 


x 


dx — 7T 2 — Z7T log( — 1 — Xl) . 


(In the second calculation, any integrals from 0 to x\ should be understood 
as including the small semicircle below —1.) Collecting the pieces together 
and dividing by 2, the final result for L+(Ci) is 

L + (Ci) = -£(1) - £(1 + 1/xi) + 7T 2 + nr log [(l+xi)/xi] . 

The calculation of £+(£ 2 ) is similar, save for the switch to Log C2 and the 
fact that the semicircular part of the contour must now go above the point 
x = — 1 instead of below it, causing + 7 r 2 to be replaced by — 7 r 2 . Adding 
everything together, 

2 2 

L +( C °) = - + £(l + l/a+)) + - 5 - log [(l + xi)(l + x 2 )x 2 /xi] . 


a= 1 


a =1 


The stationary Y-system, satisfied by the l/x a , can be used to show that the 
argument of the logarithm is equal to 1 , and hence the imaginary part of the 
sum vanishes. For the rest, substituting x a = l/y£ 3 \ and using £(1) = vr 2 / 6 , 
and the sum rule ( 3.14 ), gives 

c ( 0 ) = —20/7, 


which is the expected answer. 
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3.2 Singular lines and branch points 


The basic TBA continues to apply so long as none of the singularities of 
L\ or L 2 touch the real axis. To put this consideration into a more general 
setting, imagine that the full Riemann surface covering the complex r plane 
is marked with a set of lines, along each of which a given singularity in L\ or 
L 2 has a vanishing imaginary part. As one of these ‘singular lines’ is crossed, 
singularities change their status from inactive to active or vice versa, and 
a transition is induced in the TBA equations. There is a distinction to 
be made between the two different types of singularity that the L a {9) can 
have, one stemming from z a (9) = 0, the other from Y a (0 ) = 0. Transitions 
associated with the zeroes of the z a (9) always result genuinely different TBA 
equations, and we will refer to the corresponding singular lines as being of 
‘type 1’. The same turns out not to be true of the transitions associated with 
zeroes of the Y a (9), and these singular lines will be called ‘type If’. Further 
discussion of this point will be delayed until section T4, when a concrete 
example will be examined. 

It is important to realise that these singular lines are associated with 
transitions at the level of the TBA equations, and not with singular be¬ 
haviour of the T’s as functions of r, still less with singular behaviour of 
the scaling functions. For example, if we had chosen to set up the TBA 
equations using a contour other than the real axis for all of the integrations, 
then the set of lines across which equations changed from one form to an¬ 
other would have been differently placed. Nevertheless, the singular lines do 
provide an effective way to map the full Riemann surface, and they are also 
very useful as a means to organise the various generalised TBA equations. 
However, even for a model as simple as the one under consideration here, 
it is by no means an easy task to find their pattern. The main problem is 
that at a point on the Riemann surface near to a singular line, a singularity 
in L\ or L 2 will inevitably be near to the real 9 axis, and this tends to 
destabilise the iteration schemes used to solve the equations. This makes it 
rather hard to disentangle the genuine boundary of applicability of a given 
equation from its ‘region of stability’, the latter region depending on the 
particular numerical method employed and therefore having no intrinsic in¬ 
terest. One technique to get around this difficulty is to locate the singularity 
positions with high accuracy in those regions where iteration does work well, 
and then to extrapolate these positions into the ‘difficult’ regions around the 
singular lines. Especially in situations where data can be obtained on both 
sides of a given line, this allows its location to be found with a fair degree 
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of confidence. However the method is distinctly tedious to implement, and 
there is always the possibility that a crucial singularity has been missed. In 
some cases, it is possible to do much better. The key idea is to make use of 
the shifted conjugation property ( 2.11 ). If a singularity in L a (r,9) crosses 
the real axis, the corresponding point for L a (r,9) will be crossing the line 
Im 9 = T 7 t/ 10 , a harmless piece of behaviour from the point of view of the 
equations for e a . This makes it possible to find the functions Y a at many 
points where the direct attempt at numerical solution of the relevant TBA 
fails, simply by solving the equations for Y a instead. 

This remark can be put to immediate use, to resolve an apparent paradox 
about the behaviour of the solutions discussed so far in the neighbourhood of 
the point p=po on the positive-A line. Recall that along the initial segment 
of the positive-A line, the singularities of the L a (9) nearest to the real axis 
were at 9 = ±6a \ and all had imaginary parts equal to ±37 t/20 . At first 
sight this is puzzling: in the case of the Ising model the branch points in 
the scaling functions can all be associated with the ‘pinching’ of the real 
9 axis by pairs of colliding singularities (see |>] for a discussion), and yet 
here, as we approach the point p=po at which there is a branch point in 
c(r), the singularities in L\ and L 2 show no inclination at all to pinch the 
real axis. The resolution is rather appealing: the interacting nature of this 
model as compared to Ising, reflected in the fact that a TBA equation must 
be solved before the pseudoenergies can be found, serves to ‘desingularise’ 
the pinch. Instead, in the immediate vicinity of p=po, the branching of c(r) 
is due entirely to the multivalued nature of the solutions to the basic TBA 
equation. Get a little further away, and singularities start to cross the real 
axis as the point p=po is encircled, making the analysis similar to that of a 
pinch singularity. This fact, coupled to the severe instability of the ‘excited’ 
solution near to the branch point, makes it easy to confuse the situation 
with a genuine pinch singularity. However, with the help of the conjugation 
trick just described, the numerical problems can be circumvented and the 
ambiguity resolved. We first observe that as p increases beyond po, the 
functions Y a (9) found from ground-state solution to the basic TBA cease to 
be self-conjugate. Hence the conjugated functions Y a (9) start to differ from 
the Y a (9), and provide alternative solutions to some set of TBA equations. 
For po < p < 2.8, none of the singularities in L\(9) and 1*2(9) have crossed 
the real 9 axis, and so the e a (9) continue to solve the basic TBA equation, 
rather than one of its generalisations. We used these functions in (2 A) to 
find a scaling function in the segment po < p < 2.8 of the positive-A line, 
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with results which matched perfectly with the energy level missed by the 
ground-state solution. (Note that while the two levels appear to sit on top of 
each other for p > po on figure [3], this is an illusion caused by the equality of 
their real parts - the imaginary parts are different.) Moving beyond p = 2.8, 
pairs of singularities in the conjugated functions (first in L 2 , and then, at 
around p = 4.2, in L\) cross the real axis and the basic TBA together with 
the conjugation trick provides us, for free, with solutions to a couple of 
generalised TBA equations. 

In the ultraviolet limit, that is in those regions of the Riemann surface 
covering the neighbourhood of the point r=0, another simple observation 
constrains the pattern of singular lines rather strongly. The idea is to ex¬ 
ploit the control over the ultraviolet singularity positions given to us by the 
splitting into kink systems, and the consequent formula ( |3.7| ). This informa¬ 
tion can be used in the following way. Since log(l/r) = log(l/|r|) — 2 arg(r), 
changing the modulus of r, once the kink systems have formed, does not 
affect the imaginary parts of the singularity locations. However, a change 
in the value of arg(r) does have an effect, and may result in a singularity 
hitting the real 6 axis. Once this has happened, varying |r| again, while 
holding arg(r) fixed, just slides this singularity along the real axis without 
shifting its imaginary part away from zero. Hence the asymptotic pattern 
of the singular lines near to the origin is a collection of rays emanating from 


r= 0. Their directions follow from (3.7): for the type I lines, associated with 
zeroes of z a (0), they are 


arg(r) = ± ^Imc® ± — ^ 


(3.20) 


while for the type II lines, associated with zeroes of Y a (0), the relevant 
directions are 

arg(r) = ±Imc® . (3.21) 

For the ground state pseudoenergies solving the basic TBA equation, all of 
the Ca ^ have imaginary parts equal to 77t/10. It follows that there are no 
singular lines approaching the origin in the sector — 7t/2 < arg(r) < n/2 of 
the first sheet of the Riemann surface. However for the excited solution, 
which was initially discussed only on the positive-A line arg(r) = 77t/20, we 
saw that Imc^ vanished, and so ( 3.20| ) implies that a type I singular line 
approaches the origin on the second sheet along arg(r) = n/5. As this line is 
traversed, a zero of Z 2 ( 0 ) crosses the real axis and the pseudoenergies start 
to be determined by the first generalised TBA equation, to be discussed 
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shortly. In fact we already know a kink solution to this equation: it is 
just the excited kink solution to the basic TBA, shifted yet further in the 
imaginary 8 direction as appropriate for the diminishing value of arg(r). To 
find out what happens to this solution as the ultraviolet region is left, we 
now turn to a direct analysis of the generalised TBA equation that it solves. 


3.3 One pair of active singularities 

As r crosses the type I singular line on the second sheet, the singularities 
in L^{ 8 ) at it#® cross the real axis, the imaginary part of 8 ^ becoming 
for the first time positive. Contours in integrals involving L 2 { 8 ) must be 
distorted, and returning them to the real axis induces extra residue terms 
in the equations. Just as described in ref. |5| for the SLYM, these residues 


can be found via an integration by parts. Equations (|2.3|) and (2.4) become 


Q (8 — 8^) 2 

e a {8) = m a r cosh 6> + Log c —-^-^<£ ab *L b (0) 

Sa2(0 + 8 f) ti 


(3.22) 


and 


12 r / n , 3 

c(r) = - 1 m 2 sinh 8^ H— k'Y' / d8 m a r cosh 8L a (8) (3.23) 

7T 7T J- oo 

Cl— 1 

respectively. These equations must be handled with some care. A continuous 
branch of the logarithm, signalled by Log c , is implied. We set 


Log c S ab (+oo) = 0 ; - 4 (+oo) = 0 , 


(3.24) 


and then understand Log c (5/5') to be continuous as 8 varies along the 
real axis. For the moment we shall also assume 0 < Im 8 ^ < vr/5; the 
effects of (9^ straying beyond this strip will be examined later. Using 
4>ab = -ijjo log Sab(8), we have 


/ OO 

(f>ab(0)d8 

-OO 


—2mN a b = —2m 


1 

2 



Given the branch choice made above, and also the symmetry of S a b(8), this 
implies that 


Log c S' a& (0) = mN ab ; Log c S’ afe (-oo) = 2mN ab . (3.25) 
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Substituting 0 = into (|3.22|) , and imposing £ 2 (0 2 Uj ) = ^ to ensure that 
22 ( 0 ) vanishes at that point, then gives the final ingredient: 

2 

0 = m 2 r cosh + ^*(-^22-1) - Log c 5 2 2(202°■*)-X! ^ 26 *^ b (^ 2 °^) • (3-26) 

6=1 

The term 7r*IV 22 comes from the consistent application of the branch choice 
(13.24 ) to Log c S' 22 (0). Notice that it is equal to 37rz and not the m that 
might at first sight have been expected - the importance of this will become 
apparent during the calculation of c(0) below. The remaining logarithm 
must also be taken continuously, starting from the branch found when 9^ 
crosses the axis. This is not completely trivial, as 2 9^ may ultimately leave 
the strip |Im20®| < 7r/5, making it necessary to keep track of the sense in 
which the pole in 5 22 (0) at 0 = z7t/5 is encircled. 

All of this makes the analytic treatment a little tricky. Perhaps surpris¬ 
ingly, it is much less of a problem for the numerical solution of the equations. 
We used the following method: start by choosing an initial value for 0®, 
and then iterate ( 3.2 2| ) until a solution e a (9\9^) is found. At this stage 
the first potential difficulty arises: unless otherwise instructed, Fortran will 
always choose the principal (discontinuous) branch for the logarithm dur¬ 
ing these iterations, — ir < Im log(SyS') < ir. But any 2m discontinuities 
that this induces in the resulting functions e a {9\9^) leave the values of 
the L a = log(l+ exp(—e a )) unchanged, both in ( |3.22| ) and ( 3.23 ), and so 
this problem can be ignored. The functions e a {9\9^) are then inserted into 

( 3. 26| ) ; the error in the initial 9^ will be reflected in the degree to which this 
equation fails to hold. Again, one might worry that Fortran will spoil things 
by picking the wrong branch for the logarithm; this time, we sidestepped 
the problem by exponentiating the equation. Once this has been done, an 
adeptly-chosen portion of S’ 22 (20^) can be taken over to the left-hand side, 
and then inverted to obtain an improved estimate for 0^. The whole pro¬ 
cedure is then iterated until a stable solution is found. In this way we were 
able to map out large areas of the complex r plane, but not, as it happens, 
a segment 0 < r < r c of the real axis. The reason for this difficulty will be 
explained later, but before that we will discuss the infrared limit. 

With a little practice it is not too hard to guess how the solutions will 
behave in this regime, but just to be sure we tracked a solution and its 
singularities as r varied along the line r = 2z+(l—0.5 i)t, using the basic TBA 
( pk3 ) for 0.1 < t < 1.5, and the generalisation ( |3.22| ), ( p.26| ) for 1.5 < t < 4. 


(oh 
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We found that the singularity positions are indeed smoothly continued by the 
new equations. When r finally reaches the real axis, the singularity positions 
regain the standard conjugation symmetry under 6 —> 6*, as oppose to the 
shifted conjugation symmetry under 0 —> # associated with the positive-A 
line. However the way in which the symmetry is realised is rather different 
from the situation before any continuation, when the basic TBA equation 
applied and singularity pattern was as shown in figure 0. Then, all of the 
Pa lay on a fixed line of the mapping 9 —> 0*, namely Im# = 77 t/ 10. (Recall 
that all functions in the game are 77rz/5-periodic, so this is indeed a fixed 
line.) This time, and (3 ^ are instead swapped with their negatives 
under the conjugation; they are therefore both purely imaginary. Figure || 



Figure 8: Contour plots of \z a {0)\/(l+\z a {0)\) in the complex 9 plane, for 
the first generalised TBA equation with r= 3. Labelling as in figure 4. 

shows the functions at r = 3; note how, in contrast to the earlier pictures, 
the zeroes of lj and at /?} 0 ' and (signalled by the dots • on the 

imaginary axis) are clear of the fixed lines of 9 —> 9* and 6 —> 9. 

Once these qualitative features have been understood, it is possible to ex¬ 
tract some exact asymptotics. On the real axis with r larger than r c , # 2 ^( r ) 
lies on the segment of the imaginary 6 axis between m/IQ and m/b. Fur- 
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thermore, during the continuation from the ultraviolet end of the positive-A 
line, $ 2 °^ reaches this segment from large positive values with no monodromy 
around the point nr/10. Combining this information with the form of 822 (d) 
shows that, for real r > r c , 

Log c 5 2 2(202 O) (r)) = log5 2 2(2^°' , (?’)) + 2m . 

Equation ( |3.26| ) therefore becomes 

2 

0 = m 2 r cosh^^ - log *S' 22 (2^ 0) ) - ^0 2ft *L 6 (0^) . (3.27) 

6=1 


Since 5 22 ( 26) is real and positive for imaginary 9 between nr/10 and nr/5, 
all three terms in this equation are real. As r grows, the first two terms 
must cancel, as the convolution term tends to zero. To achieve this, <4° 
approaches the limiting value of nr/10, at which 5 22 ( 29^) has a pole. The 
leading behaviour of c(r) follows immediately: substituting 9^ = nr /10 
into ( 3.23| ) and using the relation 


7r 


mi = 2 m 2 sin(—) 


(3.28) 


gives c(r) ~ —Omir/ir, or E(X,R) ~ -Eb u ik(A,-R) + M\. This confirms that 
the state under discussion is a one-particle state, with mass gap M\. Since it 
also has zero momentum (the total momentum of a level is conserved under 
analytic continuation), it must be the first line above the ground state in 
figure |||. As a further check, a couple of corrections can be worked out. 
These should reproduce the formulae for finite-volume mass shifts found in 
18], and summarised in appendix [A]. The more precise cancellation 


refs.] 


of the first two terms in (3.27|) requires 


& 2 °\r) ~ 1 [tt /10 + tan(^) tan( 


2vr 


77127* COs(7r/l0) 


(3.29) 


This improved estimate for 9^° gives one correction, the so-called fr- term, 
when substituted into ( 3.23| ). A second comes on expanding L a (9) in ( 3.23| ), 
using the leading behaviour of e a (9) obtained from the first two terms on 
the RHS of ( |3.22|) . For this correction we can set 9^° = nr/10, and then use 
the identity 

+ =Sol ((,+ ') . (3.30) 


s.2(e - !*) 
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Putting these pieces together, the two leading IR corrections are 


(r) ~ —— (mi + 2 m ,2 cos(-^-) tan(^) tan(^-) 2 e m 2»-cot(7r/io) 


7r \ 


10 

roo 


10 


I /* OO n /JT- v 

- — Xj J dd ™a cosh 9 Sia(8 + y)e-™“ cosh0 ) . (3.31) 


These agree with the expected results, as given in ref. [18]. 

In the opposite, ultraviolet limit the singularity pattern is rather differ¬ 
ent. As remarked earlier, once the limiting form of the excited solution to 
the basic TBA on the positive-A line near to r = 0 is known, the fact that the 
only r-dependence of the kink solutions is in the ± log(l/r) anchoring points 
for the left and right kink systems allows the solution to be continued back 
to the real axis, thereby obtaining a kink solution to the generalised TBA 
equations ( 3.22 ) and ([3.26 ). Recall from just after (3/7) that for the excited 
solution to the basic TBA all of the ci ' 1 had imaginary parts equal to 77 t /10 


apart from c^\ 


which had an imaginary part equal to zero. Substituting 


into (3/7) shows that for real r the kink limit of this solution to the gener¬ 
alised TBA equation ( 3.22Q has the imaginary parts of all of the [3a\r ) equal 


to (±)77t/10, apart from (3^\r), for which Im(r) = 0. This also holds 
before the kink limit is taken, since when r is real, the singularity pattern 
must be symmetrical about the real 0-axis. For a zero of Y to move away 
from a fixed line of 9 —> 9*, there would have to be a corresponding zero, 
with equal real part, moving in the opposite direction. This is only possible 
once the left and right kink systems have met. (Numerically we found that 
this happens at the previously-mentioned point r = r c ~ 2.6646510318(2) .) 
Figure || shows the solution at r= 1. Notice that the pattern of singularities 
is indeed symmetrical under 9 —* 9* , but in a different way from the infrared 
solution shown in figure ||. The situation is similar to that observed for the 
first excited state of the SLYM in refs. ||, as r increases from the ultravi¬ 
olet, the two kink systems come into contact and the TBA profits from the 
extra possibilities that this offers for satisfying the conjugation symmetry, 
and rearranges the singularity pattern into a form more appropriate for the 
infrared limit. For real r < r c we have Re (3^ ^ 0 and Im / 3 ^ = 0, while for 
real r > r c . Re/?^ = 0 and Im E (—7r/5, 0). (There is also a transition 
associated with (3^ , when it moves from the line Im/3^° = 77t/10 onto the 
imaginary axis. This happens at r ~ 1.5, but since the associated singular¬ 
ities are inactive, the equations do not change.) The transition at r=r c will 


25 













Re 

0 12 3 4 



Re 

0 12 3 4 



Figure 9: Contour plots of \z a [0)\/{\+\z a (0)\) in the complex 9 plane, for 
the generalised TBA equation with r= 1. Labelling as in figure 4. 


be analysed further below, but first we will discuss the ultraviolet behaviour 
of the equations, where the singularities follow the pattern predicted by the 
excited solution to the basic TBA equation. 

There is a new complication in the ultraviolet regime when r is real. 


Since (ffl is real, (3.2) implies that singularities at i n -^l and L 2 

are sitting exactly on the integration contour. For analytic work, the most 
convenient approach seems to be to take these singularities into account 
via small distortions of all contours of integration, below the real axis near 
to 6 = — and above the real axis near to 9 = (3^. This is rather 
like the usual he’ prescription in quantum field theory, and amounts to 
giving r an infinitesimal positive-imaginary part. Here, the justification for 
the prescription comes from the analytic continuation idea that has been 
fundamental from the outset. (In ref. J5|, we instead rewrote the equations 
using principal values for all integrals. Whilst this puts the equations into 
an appealingly symmetrical form, the branches of the logarithms are hard 
to disentangle, and the treatment below seems to be more promising for 
systematic work.) 

To verify that the value of c(0) given by ( [3.23 ) and the solution to the 


,(0) 
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generalised TBA equation is the same as came via the excited solution to 
the basic TBA, we start by finding the kink form of the generalised TBA 
system. From (3.7) we have, as r —> 0, 

4 0) (r) ~log(l/r) -4 0) 

where, for the solution under discussion, c|4 is a real constant. Hence 

4° } ~ log(l/r) - 4 0) + y . (3.32) 

□( 0 ) - 


The le prescription for r means that 9\ is still (just) inside the strip 0 < 
Im 9^ < 7r/5. We can take this into account by setting r = pe te , with p 
real and e understood to be vanishingly small and positive. The kink form 
of equations ( |3.22| ), (|3.2(j) and ( 3.23 ) follows on shifting 9 to 0+log(l/p), 
defining £ kmk (0) = £ a (0+log(l/p)), employing the limiting form of 9^\r), 
and finally taking the p —> 0 limit, dropping inessential factors of ie : 

2 

£ a mk (#) = i mare e + Log c S , a2 (6i+c (0 )_f+?,e) - Y 4>ab*Lb mk (9 ); (3.33) 

6=1 

2 

o = im 2 e- c(0 W 5 + m(N 22 -l) - Y </, 2b *L kink (- C (°)+f-*e); (3.34) 

6=1 

/? /? 2 n qq 

c(0) = — im 2 e c< '° ) e m / 5 H— 2 Y1 / d9 ^m a e e L kmk (0). (3.35) 

(For the last equation, use is also made of the fact that half of the integral 
contributing to c(0) comes from the left-hand kink system, and half from the 
right.) Differentiating ( [1.33 ) with respect to 6 and substituting the resulting 
expression for \m a e 6 into ( 3.35|) gives 


c(0) = —im 2 e~ e 

IT 


in/b 




r d 


,_J0 


a= 1 


-*2(»+4“ ) -f)+5:Ki,*h" ,k (#) 


L 


kink 


(9) 


6=1 


The use of ( 3.34 ) simplifies this to 


’(0) = 12(iV 22 -l) + -4 Y [ Lo §£ a ( 1+e £ ) ds 

71 a=l J£a 
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7T 


E 

a, 6=1' 




<5M ink w 


L 


kink 


(*) 


= 12(N 22 -1) + 


E / LogfiaCi+e £ )* 

_ 1 J Sn. 




~ ~ 2 E \L kink (9)L ki * k (9) 


OO 

—oo 


a,6=1 

Here £ a is the contour swept out by e kmk , and Log £a is a continuous choice 
of the logarithm along this contour. The branch choice made previously 
implies that L| mk (+oo) = 0, while for L kmk (—oo) the 6 —> —oo limit of 

( H )> 


e kink (-oo) = 2mN a2 + £ iV ab L kink (-oo), 

6=1 

can be used. Rearranging and changing variables in the integrals, 

2 


(3.36) 


c( 0 ) = 12(iV 22 -l) + E^ [L + (C a )-TiiN a2 L kink (-oo) 


< 2=1 


6 z 

-12-e| ink ( -oo) H—j E L+(C a ), 

7T 7T Z J 

( 2=1 


(3.37) 


with L + (C a ) as dehned in ( |3.9| ). The contours are a little different from the 
previous case: a single zero in z 2 {9) crossed the real axis in the continuation 
from the basic TBA to the new equation. With no further numerical work 
we can deduce that this time both contours go below the point x = —1. As 
a result, ( [3. 16 ) is replaced by 


Log Cl (l+xi) = log 11 T X \| nr ; 
Log C2 (l+£ 2 ) = log |1 + x 2 \ - nr. 


(3.38) 


The branches of Log Ca (x) should be unchanged, and indeed substituting 
(3.38) into the imaginary part of ( [3.36 ) shows that (3.19) still holds. From 
this point on the evaluation of the integrals runs essentially as before, though 
note that this time the residue terms do not cancel but rather reinforce each 
other. The result is 


E L + (C a ) = - ^(£(l)+£(l+l/x a ))+7r 2 + ^log [(l+a7i)(l+x 2 )/(xix 2 )] 


a= 1 


a=l 
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Substituting into ( 3.37| ), and recalling that £ 2 mk ( —00 ) = — Log Cz (x 2 ) = 
— log(— X 2 ) + ^7^, the imaginary parts again cancel and we finally recover 
c(0) = —20/7, just as before. Assuming that the analytic continuations had 
been performed correctly, this calculation had no option but to give the right 
answer. We included it anyway to show how the generalised equations can 
be tackled directly, so long as the branches of the logarithms are handled in 
a consistent way. 

From a numerical point of view, there is a difficulty in the ultraviolet 
regime. The zeroes in the Y’s at T/?® entail divergences in the pseudoen¬ 
ergies on the real 6 axis, along which the TBA equation is being solved. As 
a result, direct numerical iteration of the equations, as outlined above, be¬ 
comes unstable as r approaches the segment [0, r c ] of the real axis. Taking a 
hint from Q, a useful trick is to remove these divergences explicitly, before 
proceeding to a solution of the equations. It turns out that this can be done 
in a particularly nice way. The first step is the following identity, a simple 
consequence of equation (9) of ref. jl^] : 


M&) = -MO) +X>ab* M&) 


(3.39) 


6=1 


where 4>h(0) = h /(2 cosh 4^) and h is here equal to 5. (For a discussion of 
the Tv-related cases of such identities, see ref. [13].) This can be integrated, 
using 


o ^ fiQ 

MO) = 0n%—Log c a h (e ± ^) , cr h (9) = tanh(—) , (3.40) 

to yield 

2 

±Log c S ad (6) = Log c a h (9±^)l [ ^ ] -J2M*Log c v h (0±^)l [ M ] . (3.41) 


Substituting the d =2 case into ( |3.22|) , recalling that (3^ = 6^—nr/h, and 
absorbing the extra terms by defining 

MO) = MO) - Log c (cr/j(6»-4°V/i(^+4 0) )) ( 3 - 42 ) 

and 

L a (0) = Log c ( a h M(3i O) )a h (0+(3i O) ) + e~^ e) ) , (3.43) 
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the excited-state TBA equation has a form identical to the original equation 
for the ground state: 


2 

£ a {0) = m a r cosh 6-^2 ct>ab*L b (9) . 

6=1 


(3.44) 


The different singularity structure is now entirely encoded in the definitions 
of e a and L a , which, in contrast to e a and L a , are not themselves singular 
at i/?®- The expression for the central charge also has a formal similarity 
with the original TBA: the integral 


J d9e lkd log (a h (0-x+t^)a h (ff+x-t^)^ = 


i sin kx 


k cosh(7r&:/h) 


(3.45) 


(convergent in the region |Irn k\ < h/2), taken at k = i and x = 9^° , allows 
equation (ft. 23) to be rewritten as 


c(r) 


— / d6m a r cosh 6L a {6) . 


(3.46) 


These equations look just the same as the original TBA system. However it 
is important that they contain an extra parameter, namely the value of . 
When handling the system numerically, we can’t search for the singularity 
at - it has been explicitly removed - but must instead focus on 9^\ 
and impose 

e ?2 ( 0 2 O) ) = zcot h(|4 0) + f) e £2(e 2 0)) = —z CO th(f^o) + f) . (3.47) 

This method was used to obtain table [I], and also figure [9]. 


3.4 Type II singular lines 

The problems encountered at the end of the last subsection can be traced 
to the presence of a type II singular line, along which a zero of Y a {9) has 
vanishing real part. As hinted earlier, the behaviour of the generalised TBA 
as such a line is crossed has some novel features. 

The first point to note is that any initial crossing of the real 9 axis by 
a zero of Y\ or I 2 must be preceded by the crossing of a zero of either 
z\ or Z 2 ~ this follows from (|3.2| ) and ( |3.3[) , and the configuration of the 
singularities before the continuation begins. Therefore, just prior to the first 
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crossing of the zero of Y\ or Y2, there will already be an active singularity 
in the equations, a zero in z c (9), say, at Without any significant loss of 
generality, suppose that this singularity became active as it crossed the real 
9 axis from below, so that it currently lies in the strip 0 < Irn O^P < 7t/ 5. 
The accompanying zero or zeroes of the Y a , at (3c J = 9^ —iir/5, lie below 
the real 9 axis and at this stage are inactive. The generalised TBA equation 
for e a {9) contains a term 

Log c S ac (9 - 9&) (3.48) 

which can be traced to the presence of the active singularity at 9 C . There 
is also an additive contribution to c(r), equal to 

— irricsinhd^ . (3.49) 

7r 


Now consider what can happen to these terms as r varies. The simplest 
option for change is for Im 9^ to become once more negative, meaning that 
a type I singular line has been crossed. Clearly this corresponds to undoing 
the trapped piece of contour, and the net effect is for the two terms ( |3.48 ) 
and ( 3.4S| ) to drop out of the equations, returning them to their original 
state. Note that this change in the generalised TBA equations is signalled 
by /3c ^ leaving the strip |Im < 7r/5. 

More subtle is the behaviour when Im 9^ grows beyond 7 t/5, pi' cross¬ 
ing the real axis but remaining inside the strip |Im (3\ < 7r/5. This corre¬ 
sponds to the crossing of a type II singular line. (We already discussed a case 
with Im 9^ exactly equal to 7r/5, found when r hits the segment 0 < r < r c 
of the real axis from above. The question currently being addressed would 
arise if r was then continued beyond this segment into the lower half plane.) 


As Im 9^ passes 7 t/ 5, the singularities in the L a (9) at (3c ’ cross the real 
axis. Further terms must be added to the equations, which are conveniently 
summarised using the incidence matrix of the T 2 diagram, l . The terms 
( |3.48| ) and (|3.49D become 


,00 


and 


Log C s ac (9 - 0&) - ]T l™Log c S ab (9 - (<#>_*tt/5)) 


6=1 


— im c sinh 9^ 
7r c 


— l^rrib sinh(0L- ) —77 t/5) 
71 b =1 


(3.50) 


(3.51) 


31 











respectively. (The minus signs are there because the new singularities are 
zeroes of Yj,($), rather than z c (9).) The situation seems to have become 
more complicated, but is helped by the following pair of identities: 

2 

Sac{0 - l-K/h)S ac {6 + w/5) = II S ' ab^) 1 ^ ! 

6=1 


m c [sinh(0— nr/5) + sinh( 0 +z 7 r/ 5 )] = ^ mb sinh 9 . 


6=1 


Dehning 9c * = 9^ — 2z7t/ 5, the terms ( |3.5C| ) and ( |3 .5 1| ) simplify to 


-Log c S ac (0-9P) 


(3.52) 


and 


- im c sinh 9^ ^ 

7r 


(3.53) 


These are exactly the extra terms that would have been found had a singu¬ 
larity in L c (9), a zero of z c (9), crossed the real 9 axis, though from above 
rather than from below. From Q3.2|) or (3.3), there is indeed a zero of z c (9) 


(i 1 ) 

at 9c J , so this reinterpretation of the equation is consistent with the singu¬ 
larity pattern actually found. 

It is worth pausing to consider what has happened. The situation under 
discussion was the crossing of a type II singular line, and we have shown that 
it does not lead to a fundamentally different equation, but just a relocation 
of an already-active singularity, back into the strip |Im 9\ < 7 t/ 5 . This is in 
contrast to the crossing of a type I singular line, which always either adds or 
subtracts a singularity from the active list. It can now be deduced that any 
equation obtained from the ground-state TBA system by analytic continu¬ 
ation in r can be rewritten in such a way that the only active singularities 
correspond to zeroes of zi or Z 2 in the strip |Im0| < 7 r/ 5 . Furthermore, in 
the rewritten equations a zero of z c {9) in this strip will be active if and only 
if Pc \ the position of the accompanying zero of Yf( 0 ) and/or ^ 2(^)5 satisfies 
|Im/3c J )| < 7 t/ 5 . This observation provides the bridge between our approach 
and the methods advocated in refs. [17, |4|. 

Finally, we should examine the behaviour of the equations for the first 
excited state near to r=r c . Recall that in the far infrared, 9^\r) ap¬ 
proached Z 7 r /10 from above. As r decreases, 9^ grows and finally, at r=r c , 
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it arrives at mjh. This is necessary for there to be a continuous transi¬ 
tion to the singularity pattern observed when r is less than r c . There are 
also inactive singularities associated with zeroes of Y\ and Y 2 , located at 
±/? 2 °^ = ±(02°) — ^*7t/ 5), and as r approaches r c . these singularities start to 
pinch the real axis. This poses a mild puzzle, the exact converse of that 
discussed in section 3.2: there, a branch point of c(r) was present and yet 


there was no evidence of an integration contour being pinched, while here, 
contours are being pinched and yet (because r is real) we do not expect 
there to be a branch point. The absence of this branch point can be under¬ 
stood if we use the standard technique for analysing a pinch singularity, as 
described, for example, in chapter 2 of |Hj. This involves continuing r in a 
complete circle around the potentially singular point, and then comparing 
the result with the uncontinued situation in order to find the discontinuity 
over the putative branch cut. 

Near to r c , we found the following dependence of O^/' on r: 


@ 2 °\r) = i^r + B(r c — r) 1 / 2 + O(ro— r) (3.54) 

o 

with B = 0.3763188677(2). This gives a good control over the positions of 
the singularities, allowing their changing status to be followed as r traces out 
a small circle about r c . At first sight the situation after continuation looks 
very different from that before: the active singularities at ±# 2 °^ have left 
the strip |Im$| < 7 t/ 5, their former positions now being occupied by a pair 
of inactive singularities, whilst the two singularities at T/^ 0 have swapped 
over and become active. However, while the assignments of singularities as 
active or inactive have changed, their overall pattern has not. Even better, 
a moment’s thought shows that what has happened is just as would be 
expected, given that the type II singular line running from 0 to r c along 
the real axis was crossed exactly once during the continuation. The earlier 
discussion of such situations applies, and if the equations are rewritten as 
described above, they return to exactly the form that they had before the 
continuation was performed. Hence the discontinuity is zero, and there 
is no branch cut after all. The functions Y a (r,9) are single-valued in the 
neighbourhood of r=r c , despite the square root in ( |3.54| ), which merely 
shows that the labelling of the singularities is not single-valued as r c is 
encircled - it does not of itself require that the Y 's should be multiple¬ 
valued. 

One particular consequence of this discussion is that the type II singu¬ 
lar line actually terminates at r=r c . This behaviour is special to the lines 
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associated with zeroes of the Y’s, and is possible because the number of 
zeroes of the Y a (9) in the strip jlrri 9\ < 7t/ 5 does not change as such lines 
are crossed. By contrast, type I lines separate regions of validity of funda¬ 
mentally different TBA equations, and so they cannot just stop at isolated 
values of r. In fact, they always seem to be anchored on some covering point 
of r=0, although we cannot absolutely rule out other closed loops elsewhere 
on the Riemann surface. 


3.5 On to the third sheet 

The third sheet can be reached by analytic continuation around either B or 
B on figure [|. The point B lies in a part of the Riemann surface described 
by the basic TBA equations (2.3) and (2.4), and continuing around B in 
an anticlockwise sense we find a third set of solutions to this system. If 
r is now returned to the origin keeping arg(r) greater than approximately 
2.1 l7r/5, no singularities cross the real 9 axis and a new kink solution to the 
basic TBA is uncovered. The value of c(0) on the third sheet can now be 
calculated much as before, using the new forms of the contours C\ and C 2 


as the necessary input into the formula (3.5). We found that C\ runs from 
the origin to a positive real endpoint x\, with no winding about x = — 1, 
and that C 2 runs from the origin to a point X 2 < —1, passing above x = — 1. 
The relevant stationary solution to the Y-system is therefore the s=5 case 


of (|3TlD : 


xi = I/ 2/1 


(5) 


x 2 = I/ 2/2 


(5) 


From the contours follow the relevant branches of the logarithms: first we 
deduce 

LogCiU+aq) = log|l + xi|; 

Log C2 (l + x 2 ) = log |1 T x 2 \ + nr , (3.55) 

and then feeding this into the imaginary part of ( |3.17 ) establishes that 
Log^(xi) = log|xi| - 2 z7t; 

Log C2 (x 2 ) = log |x 2 | - 3z7r . (3.56) 

The branch of Log Cl (aq) means that, despite x\ being a positive real number, 
L+(Ci) has an imaginary part. Extracting this and dealing with the rest via 
the change of variables t = x/(l+x), 


L + (Ci) = JZ(x 1 /(1+xi)) + nr log(l+xi) 
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For L + {C 2 ) the calculation is just as described in section 3.1, modulo the 
revised branch of Log C2 (x). The result is 


L + (C 2 ) = -£(1) - C(1 + 1/X 2 ) ~ -7T 2 + — 77T log x 2 (l+x 2 ) 


Adding the two together, the stationary Y-system can be used to show that 
the imaginary part of the total vanishes, and via the sum rule (3.15) we 
obtain 

c(0) = -68/7 

which is indeed the expected answer for the third sheet, given that at r=0 
it sees the state generated by ipn. 

During the continuation around B , the singularity at /l® remains in the 
left-hand kink system, but those at (Pp and defect to the right-hand 
system. In the notation of (377), the amounts to 


r 


r 


+ 1 
- 1 
+ 1 . 


(3.57) 


Although the singularity at pp doesn’t change its allegiance, it does move 
in a significant way, and in the new kink limit it comes to be symmetrically 
placed with the defecting —ftP- The relevant data about the asymptotic 
disposition of singularities is conveniently summarised using the constants 


cP defined by (|3.7|): 


1—1 

3 

p o 

II 

0 , 

Irn 4 0) 

T 0>o) 

Imcj = 

±-, 

10 ’ 

Irn cP 



Im 


In addition 


= +1.10933(9) 
= -1.10933(9) 


= ±- 


7ir 

To 


Rec^ 01 = Rec^ 


(3.58) 


(3.59) 


The imaginary parts of the singularity positions are important because, 


via equations ( 3.20 ) and ( 3.21 ), they determine the directions along which 
singular lines leave the point r=0 on the Riemann surface. The pattern 
here on the third sheet is much richer than previously seen: in the sector 
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7r/2 > arg(r) > 0, there are three type I lines and two type II lines. Their 
directions are: 


arg(r) 

= 2.10933(9)^ 

0 

(for z 2 (9)); 

(3.60) 

arg (r) 

= 1.10933(9)- 
5 

(for Y 2 (9)); 

(3.61) 

arg(r) 

7T 

5 

(for z\(9 )); 

(3.62) 

arg(r) 

= 0.10933(9)^ 

5 

(for z 2 {9))\ 

(3.63) 

arg(r) 

= 0 

(for Hi(0)). 

(3.64) 


(Just beyond this sector there is, as always, an infinite accumulation of 
type I lines, with arg(r) = 77r/10 — 7t/5 = 7t/2.) As r moves towards the 
real axis keeping near to the origin, with arg(r) decreasing from just below 
vr/2 to zero, the singular lines are crossed in turn, and various generalised 
TBA equations are encountered. First, the basic TBA applies; then, as the 
type I line ( 3.60| ) is passed, (3^ enters the strip |Irn fj\ < tt/5. The zeroes of 
z 2 {9) at ±# 2 ^ cross the real 9 axis and we move into territory described by 
equations (|3.22|) , (|3.23| ) and ( |3.26| ). However this time it is the singularities 
at that have become active, rather than those at ±$ 2 °\ an d so strictly 

speaking the labelling in the generalised TBA equations should be adjusted 
accordingly. This has an interesting consequence. The same region of the 
Riemann surface, namely that part of the third sheet described by equa¬ 
tions (pi), (|T23|) and (lT26|) , can alternatively be reached from the basic 
TBA by first continuing around the point A, and then looping anticlockwise 
around the point B. As A is encircled, the singularities at ±9^ become 
active, and these remain the only active singularities during the rest of the 
continuation. Their ultimate positions are those at which the singularities 
at ±# 2 ^ previously been found, when the continuation had been about 
B. Therefore, if the first continuation is combined with the reverse of the 
second, the net effect is to swap over 9^ and 9^. In the last subsection it 
was observed that a singularity can be swapped with its negative under con¬ 
tinuation around certain points on the Riemann surface; here is an example 
of a more complicated exchange, associated with a continuation around a 
non-trivial cycle on the surface. 

Passing the type II line ( p.61 ), the next singularities, zeroes of Y 2 (9) at 
cross the real axis. This does not change the set of j3a^ inside the 
strip |Im/3| < 7 t/ 5, and the only effect of the transition is to reallocate the 
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active singularities. Equations ( 3.22| ), ( |3.23| ) and ( p.26| ) continue to hold, 
but now with the active singularities at ±0^ \r) = ±(9^\r)—2ir/5). 

A genuinely new equation is only found when the type I line ( 3.62| ) is 
crossed. The singularities in L\{9) at become active, and a generalised 
TBA equation involving four active singularities results. It should by now 
be clear that the labelling of singularities is path-dependent, and so for 
notational simplicity we shall opt to label the two pairs of active singularities 
as and with 9^ and 0® both having positive imaginary parts. 

This means that 0® belongs to the right-hand kink system, and 9 2 0 ^ to the 
left-hand one. The new equation is then 


e a (9) = m a r cosh 9 + Log c 


with 


<gqi(M 0) ) S a2 (9-9^) 
Saiie+9^) S a2 (9+9 ( 2 0) ) 


~^2(t)ab*L b (9) . ( 3 . 65 ) 


6=1 


3 ^ roc 

c(r) = —-* (mi sinh0[ O) + m 2 sinh^ 0) ) H— ■= / d9m a r cosh 9L a {9) 

vr 7T a=1 7-0o 


( 3 . 66 ) 

There are now two independent singularity positions to be fixed, and the 
equations for thes 
(|3.65|) . This gives 


equations for these are found by setting 9 equal to 0, and then 9 2 in 


in = m a r cosh 9^ + E Log c 


c ( a(0)_a(0)\ ^ 

r ^ l - e ^»^ o) ). 0=1.2. 


6=1 


Sab(0, 


< 0) +9‘ 0) ) 


6=1 


( 3 . 67 ) 

There are 2mri ambiguities on the left-hand sides of these equations, both 
from the branch choices involved in solving exp(— £ a (@a°^)) = —1, and from 
the memory of the particular path of continuation implicit in the notation 
Log c . As seen in the last subsection, it is important to sort these out if a 
direct analytic treatment of the ultraviolet asymptotics is to be given. We 
will leave further discussion of this point for now, contenting ourselves with 
the calculation that was performed earlier, in the sector of the third sheet 
where the zero-singularity TBA applied. These worries are unimportant for 
numerical work, and using the exponentiation and inversion idea described 
earlier, equations ( |3.65 ) - ( 3.67 ) can be solved over large regions of the 
complex r plane. In particular, they turn out to hold even when r is real, 
so long as the infrared regime has been reached. However near to r=0 
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there is one more transition to go, caused by the type I line ( |3.63| ). As this 
singular line is passed, the zero of I 2 (0) at (3^ leaves the strip |Im/3| < 7r/5 
and the associated singularity ceases to be active. There is now only one 
pair of active singularities, and in contrast to the previous example, these 
singularities are in L\ and not L 2 . The relevant equations can be found 
from (|3.65) - (3.67) by deleting all terms involving 6^: 


e a (9) = m a r cosh 6 + Log c 


s a i(e-e[ 0) ) 

s a i(o+e[ 0) ) 


~ 0ab*-k&(0) > 


6=1 


c(r) = — 


12r z sinh 0| O) + 4 ^ 


7 r 


7T 


2 

a =1 ' 


d9m a r cosh QL a (9) , 


(3.68) 


(3.69) 


and 


m = m\r cosh $[ 0) +Log c ^^||^ - 4>ab*L b (9 { ^) . (3.70) 


5h(2^' 


6=1 


These equations should hold in the sector 0.10933(9)7r/5 > arg(r) > 0 while 
|r| remains small enough that the splitting of the pseudoenergies into left 
and right kink systems holds. However the proximity of the type II singular 


line ( 3.64 ), and the resulting presence of zeroes of Y\ (9) near to the real 6 
axis, means that a direct numerical treatment fails, at least with the methods 
that we have at present. Fortunately the equations can be recast into a more 
tractable form. As before, the key is to use the identity ( |3.41 ) to eliminate 
some of the singular behaviour in e a and L a . Here, the d= 1 case of ( 3.41 ) 
is relevant, and the revised definitions of e a and L a are 


e a (0) = e a (0) - ffLog c (a h (e-p{ 0) )a h (e+p[ 0) )) (3.71) 


and 


l n {6) = Log c ^(a h (0-/3j°V/ l (0+/5i O) ))'“ l21 + j . ( 3 . 72 ) 


Only £2 and L 2 are changed by this manoeuvre, as expected given 
The generalised TBA then falls into the more traditional form of equations 
( |3 .44 ) and ( 3.46| ), with of course the proviso that the singularities must 
be correctly placed. This equation turns out to be valid for real r out to 
4.7271(0); some numerical results are compared with TCSA data 


r = r, 


C2 


in table H of appendix 
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To understand why the equation breaks down at r = r C2 , we tracked the 
motion of the first few singularities as r varied on the segment [0 ,r C2 ]. At 
r = 0, their positions follow from the kink system. Labelling the singularities 
consistently with (|3.65 ), the right-hand kink system contains amongst other 
things an active singularity at 8*f\ with imaginary part equal to iir/5, and 
a pair of conjugately-placed inactive singularities at —8^ and d^K with 
imaginary parts equal to 0.10933(9)7 t/ 5 and —0.10933(9)vr/5 respectively. 
As r increases, the real part of — 8^ decreases, until at r = r ci = 4.691(0) 
it vanishes completely and the singularity at —8^\ part of the right-hand 
kink system, collides with one at —8^\ coming in from the left-hand system. 
At this stage the imaginary parts of these singularities are still positive. 
During the collision, the singularities undergo a 90 degree ‘scattering’ and 
for r just larger than r ci , they have turned into a pair of singularities on 
the positive imaginary axis, one heading up, and one heading down towards 
0=0. Choosing to allocate the label —8^ to the second of these, we were 
able to track its progress down the imaginary 8 axis for a short way before 
our numerical methods became unstable. Extrapolating beyond this point 
and also making use of data from larger values of r, we found that —8^ hits 
the origin when r arrives at r C2 . The natural explanation for this behaviour 
is that the type I line ( 3.63| ), which started at the origin moving away from 
the real r axis, has curved round and hit the real axis again. (It then returns 
to the origin on the complex-conjugated path in the lower half plane.) We 
confirmed this picture by tracking singularity positions for complex values 
of r, and conclude that the domain of applicability of equations ( [3.68 ) 
( 3.7C ) is restricted to a small neighbourhood of the segment [0,r C2 ] of the 
real axis. As this region is left, the singularities at ±0^°^ cross the real 0 axis 
and the more complicated equations ( 3.65 ) ( 3.70[ ) come into operation. 

The singularity movement described so far is plotted as the lines labelled 


Re (—02°'*) and Im (—0^0 on figure K), showing first the vanishing of the 
real part at r ci , and then the crossing of —0^ into the lower half plane at 
r C2 . Looking elsewhere on the figure, we see that Im0| 0 ^ remains equal to 
7r/5 as the point r C2 is passed. This means that the singularities at 
are still on the real 0 axis, and the transition to the infrared regime is not 
quite complete. For numerical work, we must continue to use the modified 
definitions ( |3.71|) and (|3.72| ) to cope with these singularities. Whilst it is 
only the singularities at ±/3j that have to be dealt with in this way, there 
are now also active singularities at ± 02 °\ and the equations are more elegant 
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Figure 10: Singularity movement for the second one-particle state: the 
dotted segments indicate the region where the numerics failed to converge, 
and were obtained by extrapolation from larger values of r. Transitions are 
at r Cl = 4.691(0), r C2 = 4.7271(0), and r C3 = 5.8786(1). 


if the associated points ±/3^ are ‘desingularised’ at the same time: setting 


[t 2 ] 

£a(0) = £ a {0) - Log c n (a{e-l3^ ) )o(0+$ ) )) Lad (3.73) 


d =1 


and 


Laid) = Log c (f[ (vW-pTMe+ffi)) 

\d =1 


,[T 2l _ ' 

ad _|_ e -Ea(0) 


(3.74) 


they again acquire the form of equations (3.44) and (3.46). This time, 
though, there are two independent singularity positions to fix: 


= -* cot h(§^ 0) + f) ; 

e e 2 (d ( 2 0> ) = _* coth + f) • 
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(3.75) 










These equations were used to obtain tables || and || of appendix [I|; the second 
set was also confirmed by the direct use of the system ( 3.65| ) - (|3.67|) . As 
can be seen from figure [l^, there is a final transition to the infrared regime, 
C3 ~ 5.8786(1), where the real part of 9^ vanishes and all active 


at r = r, 


singularities take up positions on the imaginary axis. 

As r grows further, both 9 [° and 9.p approach Z7r/10, and the identity 


7r 


m 2 = 2sin(—)(mi+?n 2 ) 


( 3 . 76 ) 


leads to the result c(r) ~ —6m 2 r/7r, or E( A, R ) ~ -Etmik(A, R) + M 2 . This is 
as expected, and just as before we can supplement the result by calculating 
two further corrections. For a = 1,2 we set 9^\r) = 77r/10 + <5 a (?’), where 
both <5i(r) and <5 2 (r) are purely imaginary for real r larger than r C3 , and tend 
to zero as r —* oo. Imposing ( 3.67| ), with the convolution terms ignored as 
being subleading, leads to the following pair of conditions: 

si 2 (<$i-<y 2 ) 


5n(z7r/5+25i)5i 2 (z7r/5+hi+(5 2 ) 

Snih-Si) 


~ e 


— m\r cos(7r/10) 


~ e 


— 77127’ COs(7r/l0) 


5i 2 (z7r/5+i5i+(5 2 )S , 22 (?.7r/5+2(5 2 ) 

(Note that with the equations exponentiated in this way the branch ambigui¬ 
ties have gone.) It will be convenient to have a notation for the most singular 
parts of the various poles controlling these equations: given a function f(9 ), 
we let lZf(9o) denote the leading coefficient of its Laurent expansion about 


f(9)=TZf(9 0 )(9-9 0 r+ J2 Ak(9-9 0 ) k ; Kf(9 0 ) ? 0. ( 3 . 77 ) 

k=p +1 

(If the pole at 9 o is simple, then lZf(9o) is just the residue.) The limiting 
behaviour of <5i and <5 2 is then determined by 
1 61+62 


Sh(v7t/5) 7T.S'i 2 (z7t/5) 


~ e 


—mir cos(7r/10) . 


61+62 


262 


~ e 


—77727’ COs(7r/l0) 


77,5'i 2 (77t/5) 775 22 (z7t/5) 

Solving, and using ( |2.1| ) for the S-matrix elements, we have 

6 2 (r) ~ _ i tan(—)e- (m2 - mi)rcos(ff/10) , 

5 

<5i(r) ~ 2ztan(—)tan(—) 2 e- mircos(7r/10) -<5 2 (r) 

10 5 
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Note that Z7r/10 is approached from above by 9^°\ and from below by 9^\ 
The leading effect that this has on the value of c(r) is through the term 
^*(mi sinh 9^+m 2 sinh 9.^) in (3.66). The result matches the contri¬ 
butions of the field-theoretic //-terms, given in appendix [A]. Notice that 
this time there are two distinct exponentials, reflecting the two processes 
2 —1 2 —2 and 2 —> 11 —> 2. A different type of correction comes from 
the thus-far neglected integral in (3.66|), using the first two terms on the 


right-hand side of 

3(0) 


to estimate L a (9). The modifications to 9 j 0 ' and 


/ 2 have a subleading effect on this term, and so their asymptotic values of 
Z7r/10 can be used, whereupon the following identity becomes relevant: 


s al (e+%) SrtV+'fi) a 

c (A ITT \ Q (A 27T \ ^a2\y \ ~ ) 

bal\V iq / 10 / ^ 


(3.78) 


The final result is 

c(r) ~ —(m 2 + 2C 1 m 1 e- mirsin ( 2n W + C 2 m 1 e- rnirsin W 5) 
7r V 


1 

27T 


^ r d9m a cosh 9 S 2a (0 + —)e- rmaCOshe ) . (3.79) 

a =1 J ~°° 2 ' 


where 

7T 3tT 27T 

C\ = 2 cos(—) tan(—) tan(— ) 2 = 24.79837... 

10 10 5 


(3.80) 

and 

C 2 = —2 sin(—) tan(—) = —3.61803 ... 

5 5 

(3.81) 

and the identity 


(m 2 — mi) cos(7r/10) = mi sin(7r/5) 

(3.82) 


was used to rewrite the second //-term. These match the predictions from 
field theory recorded in table 2 of or appendix [A|. 


4 Generalisations 

In this section we indicate how some of the above results for the T 2 model 
can be extended to the rest of the TV series. Most of our analysis is confined 
to the infrared region, where an elegant pattern emerges. Even for real r, 
the story becomes more complicated outside this region, as the movement 
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of singularities causes the equations to change their forms. The simplest 
examples were seen in refs. ||; || and also in section T3 above, with the 
single transition that could be traced to the arrival of a pair of singularities 
exactly on the real 9 axis. As was seen in that section, the effect is rather 
mild - in particular, the infrared equations continue to apply if r is given 
an infinitesimal imaginary part. The case studied in section |3.5| , the second 
1 -particle state of the T 2 theory, provides the first example of the more 
complicated behaviour that we expect to see in the general case: one pair 
of active singularities, present in the infrared, disappears completely in the 
ultraviolet. For general N we will see that the inverse process is also possible: 
a certain number of new singularities can become active at small r. This 
will be described at the end of the section, using an example from the T 4 
theory. 

Our initial conjecture is based on the behaviour already observed for 
N =1 and N=2, and the following two identities. If we set 


B = N — A + 1 


then 


and 


N 


7T ^—> 

2 sin ^ mb = mA 


N 

n 

b=B 


Sab(9 + l 


b=B 


2 h> 


7T s 


Sab(6 - 1'^h) 


— Sa<i{9 + * 2 ') 


(4.1) 


(4.2) 


These generalise equations (|3.28|) , (|3.30| ) and ( 3.76p , (|3.78 ) respectively, and 
make natural the supposition that when r is sufficiently large, the A th one- 
particle state is exactly described by a generalised TBA equation with A 
pairs of active singularities, these singularities being found in the pseudoen¬ 
ergies £b(9), b = B ... N, with all of them tending to nr/2h in the infrared 
limit. Assuming that this singularity pattern can be achieved through some 
path of continuation from the ground state, the generalised TBA equation 
must then be 


N 


£a{0) = m a r cosh 6 + ^ Log c 


(Oh 


b=B 


Sab(9-e ( b 

Sab(9+9i 0) ) 


N 


~ X] ^a6*^b( 0 ) 


(4.3) 


6=1 


with the active singularity positions 6 ^, a = B ... N, determined by the 
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equations 


m = m n r 


N Q (d 0 )o(°)\ N 

cosh e<») + £ LofeM^A) - £^.L t (#£») , (4.4) 


b=B 


s a 6 (^ 0) +^ 0) ' 


6=1 


and c(r) given by the formula 

i o r A , n , 3 A roc 

c(r) =- i(^2 m b sinhflft ) 4 —^ m a r cosh 9L a (9) . (4.5) 

71 b=B 77 a=l —o° 

At large r the dominant term is found by setting 6 = iir/2h and using 
flOl) , giving c(r) ~ —6m,Ar/ir or E(X,R) ~ ^buik(A, R) + Ma■ This is 
correct for a one-particle state with mass via- A couple of further corrections 
can be obtained exactly. The simplest to find is the A-term, and this comes 
on using the asymptotic £ a {0) ~ m a r cosh 0 — 'Log c S a A{0+nr /2) in (|4.5| ) 
(equation ( |4.2|) ha ving been used in order to simplify the ratio of S-matrices 
appearing in (|4.3| ) ). This gives the following correction to c(oo): 

3r , A r°° nr 

d9m b co S hd S A b(0 + -)e- mbrcoshd . (4.6) 

77 6=1 J ~°° 

With the leading asymptotic of the ground state subtracted off, this matches 
the F-term (A.l). The derivation of the second correction is more tedious. 
We set f9a°' 1 = m/2h + 5 a , with 5 a purely imaginary, and start from an 
infrared form of equation (|4.4|) , obtained by dropping the convolution term 
and then exponentiating: 


N 

n 


Sab(5 a Sb) 


~ —e 


—rm a cosh 8, 


( 0 ) 


a = B, B+l,... N. (4.7) 


b=B Sab(iir/h+5 a +5 b ) 

A simple consideration of the signs of the terms in these equations, using 


Sign S ab {ix) = 


(- 
(—l )' 5 ® 6 

[Tv] 


— ir/h < x < 7 x/h 
7 r /h < x < 2ix/h 


for x real, and the explicit form of l l ab , establishes that 


Sign — = (-1) 

l 


B-\-a 


a = B, B+l,... N, 
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and 


|<5b| > |<5b+i| > • • • > |<5jv| • 

However more precision is needed if we are to capture the asymptotic. Not¬ 
ing that only the S-matrix elements with ^ 0 are singular at 9 = nr/h, 
the dominant behaviour of ([O]) is 


n to + fc) 1 ™ 

b=B 



with 1Z as defined in equation (|3.77|) . Now ([h^) can be used at 6 = nr/2h 
to simplify the right-hand side, giving 


N 


n + 0 * 


I t n 1 
b ~ —e 


—m a r cos 


b=B 

Explicitly this reads 

{5 b + <5b+i) 
(5b+i + 5 b) (5b+i + 5b+ 2 ) 


W 2h )nS Aa (i!$±ir ), a = B,B+l,...N. 


_ e ~m B r cos(n/2h) n s AB (l^Tr) 
_ e -m fl+1 rcosW2fe) 7l ^^ +i ( l W ff ) 


+ 5n-2)(5n~i + 5n ) ~ (-/ 2 fc )^, iV _ 1 (*^ i 7 r ) 

+ ~ ^ e -m N rcos(n/2h) ftS A)iV (*^7r) 


Now a chain of substitutions gives 


5b + 5\ 


6+1 

5n 


\K N e-P” r , 


— B, B-\- 1,... iV—1; 


with 


and 


Vb = cos(7t/2 h) ^ (—l) 6+c m c 
c=B 


K b = n ( n (- 5 Ac (z^ 7 T 




(- 1 ) 


b+c ' 
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These values will modify c(r), at leading order, through a correction term 
equal to 


12 r 


N 


A M c(r) = - 1 E m b cos(7r/2/i)<5fc 

^ =B 
■N-i 

E + hn&n 

,b=B 
J-l 

E 2 ^ b K b e~^ r + fi N K N e~^ r 


7T 

12r 

7T 

6r 




(4.8) 


The calculation now splits into two. First, consider a term in the final sum 
for which b—B is odd, equal to 2fc—1. Then we can use the following mass 
relation: 

B+2k—l 

cos(7t/2 h) E (-l) s+c+1 m c = m fc sin(^^7r), 
c=B 

or 

MB+2/c-l = /^A,k,A-k , 

where the notation HA,k,A-k is as used in equations ( A.2 ) and (|A.5| ). Re¬ 
calling that A E{R) = —^Ac(r) = — ^^-Ac(r) and comparing with the 
formulae of appendix |A|, we see that the correction due to the a-process 
A —> k, A—k —► A will be reproduced, so long as the following identity 
between S-matrix residues holds: 


/B+2k-l 

n ( fl ( 5 A c (*^vr)) 


(- 1 ) 


S+c+1 


= ns Ak (i^n) 


\ c=B 


We checked this algebraically up to N= 5, and then numerically to A=15. 
Up to a sign ambiguity encountered when taking a square root, the identity 
can also be deduced from the following curious property of the S-matrix 
elements: 


B+2k-l 


(Sac(& + *t^-7t)SUc( 0 + *^T 7r )) 




c=B 


S A k(0 + 

SAfc(^^TT) 


If b—B = 2k is even, the story is much the same. The relevant mass 
relation is 


B+2k 


cos(7t/2 h) E (—1 ) B+C nic = mjv-fcsin(A±A_A 7r ) 


c=B 
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or 


M B+2k ~ HA,N-k,N+l-A+k ■ 

This has a chance to match the f3 correction (|A. 6 | ), from the process A 
N—k, N+l—A+k —> A. The necessary residue identity is 


/ B+2k 

k n (^ (*w 

V c=B 


7r 


(- 1 ) 


B-\-c 


= 'R-S A)N _ k (i A+ ^ k ir) 


which follows, again up to a sign, from another property of the S-matrix 
elements: 


B+2k 


II (s Ac (0 + *W 7r )' S ' Ac(9 + *%r 7r )) 


(-l) s+c S Ak {0 + 


c=B 


s Ak (e - % 


h 

A+N-k 

h 


7T) 


To finish, note that the term in ( |4.S| ) proportional to /vv receives special 
treatment, in that it is without the factor of 2 that multiplies all of the 
others. If A=2k, then b—B is odd for this term, and it corresponds to the 
a-process A k,k —► A, whilst if A = 2k+l, b—B is even and the term 
matches the /3-process A —> N—k, N—k —* A. In both cases the missing 2 
is accounted for by a symmetry factor. 

The intricate mechanism by which all of the //-terms are successfully 
recovered leaves little doubt that equations (13) • • ( fO|) describe exactly the 
zero-momentum one-particle states at suitably-large values of r. Transitions 
are inevitable as the ultraviolet regime is approached, and we have already 
seen examples of the sort of phenomena that can occur. This section ends 
with one more, showing yet another way in which the equations can change 
their form in the crossover region. 

Consider the first one-particle state in the Z 4 theory. As just explained, 
at large r we expect this to be described by a generalised TBA equation with 
one pair of active singularities, zeroes of z^O). These are at ±0 ^, with 9^ 
purely imaginary, tending to Z7r/18 from above as r grows. Conversely, as 
r decreases our numerical solutions indicate that Im Q i f ' 1 increases towards 
7 t/ 9 , until at r c ~ 2.266315266(7) there is a square root singularity in its 
value as a function of r. This transition is shown in figure [ll], and is of 
the simplest kind seen in previous sections. Decreasing r further, 
starts to run along the line Im0 = n/9, with accompanying zeroes in Y 3 and 
in Y4, both at /3^°\ running along the real 6 axis. But more importantly, 
the numerical solution reveals that a so-far inactive singularity also comes 
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Figure 11: Singularity movement for the first one-particle state in T 4 (Re 
is equal to 0 for all r in the range plotted, and so this line has been omitted). 

into play. This is a zero of zs(9), at 9^ say, which for large r approaches 
Z7r/18 from above. As r becomes smaller, Im#® at first grows, reaching a 
maximum at r ~ 1.8, and then decreases, as shown on figure |ll]. At around 
r = 1, our program became unstable, but extrapolation predicts that at 
r ~ 0.9234 the singularities at ±# 3 ^ will cross the real axis and become 
active. At this stage 9^ is still on the imaginary 9 axis, and a natural 
scenario for its subsequent movement is that it ultimately reaches —tir/9, 
undergoes a transition there and finally runs along the line Irn 9 = —7r/9. 
The ultraviolet solutions for r real would then have both (3^ and (5 4 ^ sitting 
on the real axis, K*(#) having four zeroes on the real axis, at and 

and Ys(9) having two, at Via the Y-system, the zeroes of Y%{9) also 

force Y2(9) to vanish at Unfortunately, we have not yet developed 

the computer code necessary to verify this picture directly. In its absence 
we can at least make a couple of checks. First, the zeroes just described tell 
us the signs of the numbers y a = Y^ lllk (— 00 ) for the kink versions of the 
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solution under discussion, since we already know that Y kink (+oo) = +oo: 
Sign (yi,y 2 , 2 / 3 ) 2 / 4 ) = (+1,-1,-1,+1) . 


These match the s=3, A=4 case of the solutions 


sin(sa7r/(/i+2)) sin(s(a+ 2 )- 7 r/(h+ 2 )) 
^ a sin(s7r/(/i+2 )) 2 


(h = 2N+1) 


to the stationary Tjy Y-system, which generalise ( 3.11 ). The second check in 
some senses subsumes the first. Recall from the T 2 discussions that the kink 
system for the first generalised TBA equation, studied in section was 
alternatively found using the second, excited solution to the ground-state 
TBA equation on the positive-A line. For X 4 we found a similar multiplicity 
of solutions to the ground-state equation. An excited solution to the basic 
Z 4 TBA equation can be found by continuing once anticlockwise around the 
point ro = poe llm Z 36 , po ~ 1.665(9). The limiting singularity pattern of 
this solution on the positive-A line r = p e ll7n / 36 near to p = 0 turns out to 
match precisely with the predictions just made on the basis of the assumed 
transitions in the generalised equation along the real r axis. 

Thus not only do singularities become inactive as the ultraviolet is ap¬ 
proached, as seen in section [3.5| , but also previously inactive singularities 
can become active. More work will be needed to unravel the full pattern of 
these transitions, but we have at least exhibited examples of the full range 
of possibilities. 


5 Conclusions 

The theories studied in this paper provide a useful set of toy models, on 
which it has been possible to test and develop the method of analytic con¬ 
tinuation of TBA equations advocated in [||]. 

Although relations with the sine-Gordon model at /? 2 / 87 r = 2/(2IV-|-3) 
||, and with the Z 2 Ar-|_i-symmetric conformal theories perturbed by their 
first thermal operators 10 , lend the Tj y theories some additional physical 
interest in their own right, the main objective has been to find general 
patterns of wider applicability. In the following, we summarise the main 
results obtained, and indicate some open questions that remain. 

• Shifted conjugation symmetry: The analyticity of the Y^s as func¬ 
tions of the variables a + and o_, near to a + = a_ = 0 , was conjectured 
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in section 2, and as a consequence the shifted conjugation symmetries 


and (|2.11|) were deduced. These enabled the TBA equations to 
be studied in regions where the direct numerical approach failed to 
converge. In particular, we found in this way that the branch points 
connecting the various sheets of the Riemann surface are not pinch 
singularities, but rather arise from the multivalued nature of solutions 
to individual TBA equations. This contrasts with the ‘free’ behaviour 
found for the thermally-perturbed Ising model. 


• Singular lines: Type I and type II singular lines were defined in 
section [1.2| . The type I lines divide the Riemann surface into disjoint 
regions, within each of which a different generalised TBA equation 
holds sway. The type II lines are less significant: the TBA equations 
before and after continuation through these lines can be recast into 
equivalent forms. As explained in section |3.4| , this is reflected in the 
fact that the crossing of a type I line changes the number of zeroes of 
the Y a {9) in the strip |Im0| < ir/h, whilst the crossing of a type II 
line does not. Mapping the pattern of singular lines gives an insight 
into the relationships between the various generalised TBA systems, 
and a number of techniques were explored to help in this task. In 
particular, the result described at the end of section [b^, based on the 
ultraviolet splitting of pseudoenergies into kink systems, allows the 
asymptotic pattern of lines at small r to be controlled on any sheet, 
just from knowledge of the singularity locations for the relevant set 
of kink pseudoenergies. So far we have only studied this in detail for 
the T/v systems, but it is clear that the same picture will hold for any 
TBA system which splits into separated kink systems in the ultraviolet 
regime. 


• Exotic solutions to the basic TBA equation: A basic TBA 
equation, with no active singularities, has direct access to a series of 
branch points in the scaling functions. These always turn out to be 
approximately lined up along the imaginary r axis, as explained for 
the case of the SLYM in section 3 of ref. (5|. (Figure 4 of that paper 
shows two of these points for the SLYM, whilst the points A, B and 
D on figure [l] are three examples from the T 2 theory.) Continuing r 
anticlockwise around these points, we remain in the domain of validity 
of the basic TBA, and uncover some previously unsuspected solutions. 
Furthermore, by returning r to the origin along lines of constant argu- 
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ment, we found, in sections EO and 3.5, the ‘excited’ kink solutions on 


the second and third sheets of the X 2 Riemann surface, without ever 
having to solve a generalised TBA equation. It would be interesting 
to know whether all kink solutions, both in these and other models, 
can be obtained in this way. 

Analytic continuation and sum rules: Continuation in r allows 
the contours taken in the generalised dilogarithm functions L + (C a ) to 
be determined unambiguously, and in a way which should work com¬ 
pletely generally. This has the potential to provide a more ‘physical’ 
justification for the continuations discussed as purely mathematical 


properties of dilogarithmic sum rules in |16|], and deserves to be stud¬ 
ied in greater depth. 

Desingularisation of the generalised equations: By explic¬ 

itly removing the singularities of the pseudoenergies e a (8) in the strip 
—n/h < Im0 < n/h (all related to zeroes of the Y a (6) at the same 
points) we were able to recast the generalised TBA equations, includ¬ 
ing the expressions for c(r), into the shape of the basic TBA equa¬ 
tions, albeit with redefined functions e a and L a . This was very useful 
for the efficient numerical treatment of the equations, and may be of 
deeper significance. Note that the identities (3.3£) and ( 3.4C| ) on which 
the manoeuvre depends hold with trivial modifications for all of the 
ADET -related diagonal scattering theories, and so the technique is 
immediately applicable to a number of other models. 


• Field theory and the generalised TBA: In all of the cases that we 
studied, it proved possible to extract exact asymptotics from the gen¬ 
eralised TBA equations which matched the field-theoretic predictions 
found in refs. ||], |lj|. These results relied on a number of intricate iden¬ 
tities satisfied by the mass spectra and S-matrices of the TV-related 
models. An outstanding problem is to set this into a broader context, 
and as a first step it would be interesting to know the appropriate 
generalisations to the other diagonal scattering theories. 


• The crossover region: The detailed study made in sections 3 and 4 
revealed that even for real values of r the generalised TBA equations 
can change in complicated ways as one moves from the ultraviolet 
to the infrared, typically undergoing a number of distinct transitions. 
These transitions can both increase and decrease the number of active 
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singularities, and as each is passed, a new equation has to be solved. So 
far we can only appeal to a case-by-case analysis, and more numerical 
work will probably be needed before a general picture emerges. There 
is an intriguing similarity between the movements of singularities in 
the generalised TBA equations and the movements of Bethe ansatz 


roots near to the scaling limit of lattice models, as studied in ref. |20f| ; 
it remains to be seen whether any significance should be attached to 
this observation. 

One item of unfinished business is a study of more general states, even in the 
T/v models. Based on the experience with the SLYM discussed in ref. j|], 
it is not hard to guess how to include one-particle states with non-zero 
spin: the key is to give up the 9 —> —9 symmetry of the pseudoenergies, 
and to replace each pair {9a \ —9 a' 1 } of active singularities in the equations 


previously discussed with a pair {9a\9^}, imposing 9^ = (9^)* whenever 
r is real. Note that these pairs of singularities belong to either the left or 
right kink systems from the outset, and so we expect there to be less need for 
them to undergo transitions as the infrared region is left and the ultraviolet 
approached. Preliminary studies of the T 2 model show good agreement 
with results from truncated conformal space, with, so far, no sign of any 
transitions at all, through the whole range of r. However the full story, both 
for these states and for states with higher numbers of particles, remains to be 
uncovered. Numerical problems are rather acute when trying to track more 
than one or two active singularities, and indeed it is becoming increasingly 
necessary to develop some alternatives to the rather crude iteration schemes 
that we have been using to date. 

Finally we should mention a hope for the longer term, which is to develop 
an overall understanding of the Riemann surfaces so far only tentatively 
explored. There are many structures hidden in these surfaces - in particular 
we would emphasise their monodromy groups [21] - and it is natural to 
suppose that these will be related to some of the other pieces of algebraic 
machinery that enter into the study of integrable quantum field theories. 
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A Field-theoretic predictions for infrared asymp¬ 
totics 

Finite-size corrections to the energies of one-particle states were studied 
by Liischer in ref. [4]. To all orders in perturbation theory, he found that 
the leading corrections could be expressed in terms of universal quantities 
appearing in the infinite-volume two-particle S-matrix elements. In this ap¬ 
pendix we summarise the results of Jlj] as applied to the T/v models, drawing 
also on the convenient and somewhat generalised discussion provided by 


ref. p8 |. 

If a zero-momentum one-particle state A has a mass equal to M A on the 
infinite line, then the mass of the same state on a large but finite circle of 
circumference R receives corrections from two main sources. Those from the 
first, called E-terms in [^], give a total contribution 

A f M a = -V' — V d0e- MbRcoshe (S Ab (0+nr/ 2) - 1) . (A.l) 

k 27T J_oo 

The prime indicates that terms smaller than the overall error can be dropped 
(see ref. |l£| for more discussion of this point) and the V signals that the 
principal part of the integral should be taken. (It is anyway irrelevant for the 
Tjy theories.) These terms can be traced to the interaction of the particle 
A with a virtual particle b travelling around the space-time cylinder, as 
illustrated on the left of figure O. 


Contributions of the second type (called /j-terms in []18|] ) are due to the 
processes illustrated on the right of figure [l2|, in which A splits into a pair of 
virtual particles b and c, which travel around the world before recombining 
to form A again. The result is proportional to the squares of the on-shell 
three-point couplings f A b c , and these can be determined from the residues 
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A — 

An -F-terrn 


A- 

A /x-term 


Figure 12: The two classes of diagram responsible for the leading contribu¬ 
tions to AM a ■ 


IZSAb('iUAb) °f the two-particle S-matrix elements at the appropriate fusing 
angles iU c Ah . The total shift induced by processes of this sort is 

A, l Ma = - E 7 0 ( M 1 - I M b - M c\) HAbcRAbce~^ R (A.2) 

b,c 

where 

RAbc = - lMAbcRSAb(lUAb ) , (A.3) 

and 

HAbc = Mb sin U% b = M c sin U b Ac = U . f '. U ' sin U } f c . (A.4) 

Ma 

In these definitions, Mai,c is one if fAbc / 0 and zero otherwise, and the step 
function 

( 0 if x < 0 
0(x) = < \ if x = 0 
[ 1 if x > 0 

serves to disqualify those virtual processes A —> be —> A which cannot be 
drawn as on-shell diagrams on the cylinder. 
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A-b 


(a) 



b=l,2,...[A/2] 


h-A-b 


(P) 



b=B,B+l,-N-[A/2] 


Figure 13: A set of /r-terms in the T'n theory. 


Specialising to the non-unitary Ty theories, the non-vanishing couplings 
are as follows: 

fAbe £ IR* if c = A+b < N or c = \A—b\ / 0; 
f A bc BR* if c — h — (A+b) < N ; 
f A bc = 0 otherwise, 


and the fusing angles are 


U 


\A—b\ 


Ab 




uT b) - (^'i- 


with n(A,b) = min (A+b,h—A—b). Due to the step function in (|A.2|) , only 
the contributions corresponding to the diagrams represented in figure [13] are 
non-vanishing, and these give 

■ i\A—b\ \ 

HA,b,A-b = m b smf 1 — L 7r), 

RA,b,A-b = -iRSm(i(1-^)tt), 6=1,2,... [A/2] (A.5) 

for a-processes, and 

= m b sm(^-^7r), 

= -?KS Ab (i l -^Tr), b = B,B+l,...N-[A/2] (A.6) 


l^A,b,h—A—b 

■f^A,b,h—A—b 
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for /3-processes. Note that the two types of process contribute oppositely 
to ApAf/t, since the couplings involved in the a-processes are real, while 
those in the /3-processes are purely imaginary. The total number of distinct 
contributions is equal to A, which is reflected in the fact that there are A 
independent active singularities in the relevant generalised TBA equation. 


B Numerical results 


In this appendix, numerical solutions to the excited-state TBA equations 
discussed in section 3 are compared with TCSA data. For the first one- 
particle state, table |l| was compiled using equations ( 3.44 ) ( p.47|) , and 

table |H using equations ( |3.22|) , ( |3.23| ) and (|3.26|) . For the second one-particle 
state, table |3] used using ( |3.71|) , ( |3.72|) , and tables [| and ||, ( [3.73 ) ( 3.75| ). 

Table || was also checked using (|3.65| ) - ( 3.67 ). Accuracies for the TBA 
results were estimated by varying parameters such as the discretisation of 
the 9 axis until the results stabilised at the quoted values, whilst accuracies 
for the TCSA were estimated by varying the truncation level. The first TBA 
entry in table [| was obtained by extrapolation from larger values of r, since 
our iterative routine failed to converge there. 
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D 

TBA 

TCSA 

4 0) 

0.25 

0.23816052575(2) 

0.23816052577(4) 

2.44086172191(0) 

0.50 

0.23856619676(8) 

0.2385661967(3) 

1.74875349244(1) 

0.75 

0.23957858011(5) 

0.23957857(9) 

1.34585654321(1) 

1.00 

0.24140149459(4) 

0.24140149(2) 

1.06270397320(0) 

1.25 

0.24415555107(8) 

0.2441555(4) 

0.84614229277(0) 

1.50 

0.24785826484(5) 

0.2478582(4) 

0.67206426675(5) 

1.75 

0.25241760010(3) 

0.252417(5) 

0.52676984074(1) 

2.00 

0.25764227087(8) 

0.257642(1) 

0.40064138633(5) 

2.25 

0.26326448567(0) 

0.263264(2) 

0.28454565858(0) 

2.50 

0.26896728571(2) 

0.268967(1) 

0.16239328864(9) 


Table 1: First one-particle level, r < r c . 


57 












n 

TBA 

TCSA 

■4 01 

2.75 

0.27440978264(3) 

0.274409(3) 

0.10660570253(6) 

3.00 

0.27924699310(4) 

0.279246(3) 

0.19389451747(5) 

3.25 

0.28314394898(6) 

0.283143(1) 

0.23644517085(6) 

3.50 

0.28578526162(5) 

0.28578(4) 

0.26216147499(3) 

3.75 

0.28688156615(0) 

0.28688(0) 

0.27879066486(6) 

4.00 

0.28617389043(3) 

0.28617(2) 

0.28989626682(7) 

4.25 

0.28343648042(0) 

0.28343(4) 

0.29744353109(2) 

4.50 

0.27847823059(8) 

0.27847(6) 

0.30262252477(1) 

4.75 

0.27114269310(0) 

0.27114(0) 

0.30619467150(4) 

5.00 

0.26130666042(9) 

0.26130(4) 

0.30866417603(4) 

5.25 

0.24887745081(7) 

0.24887(5) 

0.31037226542(4) 

5.50 

0.23378918672(7) 

0.23378(7) 

0.31155300241(8) 

5.75 

0.21599846946(4) 

0.21599(6) 

0.31236819753(5) 

6.00 

0.19547988038(6) 

0.19547(9) 

0.31293014812(8) 

7.00 

0.08602148615(6) 

0.08602(7) 

0.31388982426(3) 

8.00 

-0.06688310507(3) 

-0.0668(6) 

0.31410090210(8) 

10.0 

-0.50050238863(1) 

-0.500(4) 

0.31415656313(2) 

15.0 

-2.31926770945(7) 

-2.31(8) 

0.31415926412(7) 

20.0 

-5.18417167922(3) 

-5.18(2) 

0.31415926535(8) 

30.0 

-14.0517103842(1) 

-14.0(4) 

0.31415926535(8) 


Table 2: First one-particle level, r > r c . 
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m 

TBA 

TCSA 

/>?> 

0.25 

0.8095238098(2) 

0.80952380986(1) 

3.2949747170(7) 

0.50 

0.8095238271(8) 

0.8095238272(2) 

2.6018275552(3) 

0.75 

0.8095239881(9) 

0.8095239882(1) 

2.1963626206(8) 

1.00 

0.8095247257(6) 

0.8095247258(0) 

1.9086813430(7) 

1.25 

0.8095270410(0) 

0.8095270410(7) 

1.6855402842(2) 

1.50 

0.8095327731(5) 

0.809532773(2) 

1.5032248871(9) 

1.75 

0.8095448001(5) 

0.809544800(3) 

1.3490870895(8) 

2.00 

0.8095670652(3) 

0.809567065(4) 

1.2155794073(0) 

2.25 

0.8096043028(5) 

0.80960430(4) 

1.0978356190(0) 

2.50 

0.8096613348(7) 

0.80966133(5) 

0.9925341024(8) 

2.75 

0.8097417885(7) 

0.80974178(9) 

0.8973043325(8) 

3.00 

0.8098461844(7) 

0.80984618(3) 

0.8103904953(7) 

3.25 

0.8099693378(1) 

0.80996933(8) 

0.7304470394(3) 

3.50 

0.8100972214(6) 

0.81009722(2) 

0.6564058765(0) 

3.75 

0.8102034874(6) 

0.8102034(9) 

0.5873835438(9) 

4.00 

0.810246050(1) 

0.8102460(6) 

0.522609957(4) 

4.25 

0.810164185(2) 

0.810164(2) 

0.461365983(5) 

4.50 

0.809876615(7) 

0.809876(7) 

0.402917258(4) 


Table 3: Second one-particle level, r < r Cl . 
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D 

TBA 

TCSA 


-of 

4.75 

0.8092809(5) 

0.80928(1) 

0.346424(5) 

0.015804(6) 

5.00 

0.80825467(3) 

0.80825(5) 

0.29078538(8) 

0.10280614(0) 

5.25 

0.80665753(1) 

0.80665(8) 

0.23427070(8) 

0.14891716(3) 

5.50 

0.80433518(4) 

0.80433(6) 

0.17338533(0) 

0.18056955(6) 

5.75 

0.80112353(4) 

0.80112(4) 

0.09649407(3) 

0.20409098(1) 


Table 4: Second one-particle level, r ci < r < r C2 . 


60 














D 

TBA 

TCSA 


-»4 0) 

6.00 

0.7968556702(0) 

0.79685(7) 

0.0896301999(2) 

0.2223149258(0) 

6.25 

0.7913547600(8) 

0.79135(9) 

0.1501055922(9) 

0.2368163134(7) 

6.50 

0.7844612122(0) 

0.78446(8) 

0.1861508713(2) 

0.2485780614(9) 

6.75 

0.7760128576(7) 

0.7760(2) 

0.2116358291(4) 

0.2582559981(1) 

7.00 

0.7658591165(7) 

0.7658(7) 

0.2308123309(7) 

0.2663096804(3) 

7.25 

0.7538605345(2) 

0.7538(8) 

0.2457260952(8) 

0.2730722975(5) 

7.50 

0.7398899544(4) 

0.739(9) 

0.2575655998(7) 

0.2787920410(3) 

7.75 

0.7238331001(3) 

0.723(8) 

0.2670979493(9) 

0.2836579633(7) 

8.00 

0.7055886730(5) 

0.705(6) 

0.2748522116(1) 

0.2878168593(0) 

8.50 

0.6621947642(6) 

0.662(2) 

0.2864555943(6) 

0.2944542899(1) 

9.00 

0.6091397136(7) 

0.609(2) 

0.2944284302(9) 

0.2993877628(7) 

10.0 

0.4725672915(3) 

0.472(7) 

0.3039279521(5) 

0.3058454231(1) 

15.0 

-0.8450264306(8) 

-0.84(3) 

0.3136882646(6) 

0.3137048805(9) 

20.0 

-3.2169979044(2) 

-3.21(5) 

0.3141349870(6) 

0.3141351300(9) 

30.0 

-11.100815854(4) 

-1(4) 

0.3141591977(4) 

0.3141591977(4) 


Table 5: Second one-particle level, r > r C2 . 
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